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ABSTRACT 

We classify 329 late- type giants within 1 parsec of Sgr A*, using the adaptive optics integral field 
spectrometer SINFONI on the VLT. These observations represent the deepest spectroscopic data set 
so far obtained for the Galactic Center, reaching a 50% completeness threshold at the approximate 
magnitude of the helium-burning red clump (K$ ~ 15.5 mag.). Combining our spectroscopic results 
with NaCo H and K$ photometry, we construct an observed Hertzsprung-Russell diagram, which 
we quantitatively compare to theoretical distributions of various star formation histories of the inner 
Galaxy, using a \ 2 analysis. Our best-fit model corresponds to continuous star formation over the last 
12 Gyr with a top-heavy initial mass function (IMF). The similarity of this IMF to the IMF observed 
for the most recent epoch of star formation is intriguing and perhaps suggests a connection between 
recent star formation and the stars formed throughout the history of the Galactic Center. 
Subject headings: Galaxy: center — stars: RC, RGB, and AGB — stars: formation 



1. INTRODUCTION 

A growing body of evidence suggests that the stel- 
lar population of the Galactic nucleus is distinct from 
that of the bulge. Surface brightness measurements from 
the NIR to the FIR show that in contrast to the bulge, 
the central few hundred parsecs of the Galaxy are dom- 
inate d by a flat, disk-like distribution of stars, gas, and 
dust (|Keirdll99a lLaunhardt et al.ll2002t) . This difference 
is also reflected in the kinematics: OH/IR stars in the 
central hundred parsecs show higher rotational velocities 
than ex pected for an inner sin gle-component, bulge pop- 
ulation (|Lindqvist et al.lfl992f ). 

Evidence for ongoing star formation also distin- 
guishes the nucleus from the bulge. The Galac- 
tic bulge is composed primarily of an old starburst- 
like population fo r med early, some 7-10 Gyr ago 
(IZoccali et al.ll2003t Ivan Loon et alJ 120031 IZoccali et al.l 
20061 iBallero et all 12007^ 7" In contrast, the Galactic nu- 
cleus harbors substantial young and intermediate-age 
stellar populations. Intermediate-age populations have 
been most commonly inferred from broad-band photom- 
etry (iRiekd 119871; iNaravanan et aT] 119961: iDavidge et al.l 
19971 iPhilippeti all 119991 : lAlexander fc Sternberg! 119991 : 
Figer et aTTl2004f ). Star formation tracers have also 



been used to study intermediate- age populations, in- 
cluding young supergiants and lu minous AGBs, which 
trace 10 My r - 1 Gvr populations (jSellgren et al.l [19871 : 
iBlum et al.lll996bllan . and OH/IR stars, which trace 1- 
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3 Gy r populations (W ood et alJll998t ISiouwerman et all 
1999). The young stellar populations in the Galactic 
nucleus have been the subject of much recent work. 
These populations are predominantly concentrated in 
three young, massive clusters: the Arches and Quintu- 
plet clusters (Figer et al. 1999, 2002; Stolte et al. 2005, 
and references therein), approximately 30-50 pc in pro- 
jected distance from the Galactic Center, and the Central 
Cluster, located within the central parsec (Krabbe et al. 
1991, 1995; Paumard et al. 2006, and references therein). 

Star formation in the nucleus is thought to be ul- 
timately linked to the inwa rd transport of gas in- 
duced by the Galactic bar (|Morris fc Serabvnl [1996; 
iKormendy &: Kerm icutt 2004|). However, precisely how 
gas is funneled from the outer galaxy to the nucleus and 
how this process affects the resultant star formation, 
remains poorly known. At present, there is a concen- 
tration of molecular gas 2-8 pc fr om the Galctic Cen- 
ter termed the circumnuclear disk dGuesten et alJH987t 
Uackson et al.lll993l: ISerabvn et al.lll994f ). If such a struc- 
ture routinely fuels star formation in the central parsec, 
the stars formed in this region may represent an entirely 
distinct population from the Arches and Quintpulet clus- 
ters or the intermediate age populations found through- 
out the ^200 pc nucleus. Dynamical effects due to the 
supcrmassive black hole (Genzel et al. 2003, Chez et al. 
2005, and references therein) and the large stellar and 
remnant density in the central parsec may also lead to 
population differences between the Central Cluster and 
that in the larger nucleus, either prior to star formation 
or afterwards. 

A number of investigators have attempted to ad- 
dress the nature of the Central Cluster and its rela- 
tionshi p to the larger 100-3 00 pc nucleus. Most re- 
cently, iPaumarTeLaD (|2006l ) investigated the proper- 
ties of the young stellar population in the central parsec, 
spectroscopically identifying nearly 100 OB and Wolf- 
Rayet stars. They find that the majority of the young 
stars reside in two, inclined and counter-rotating disks, 
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suggesting in situ star formation in dense gas accretion 
disks. Studies of the older stellar population in the cen- 
tral parsec have been primarily limited to broadband 
photometry. They have suggested that the fraction of low 
mass stars in the population increases with distance from 
the center, perhaps due to dynamical mass segregation 
(iPhilipp et alJ 119991: iGenzel et all l2003t ISchoedel et all 
12007ft . 

A detailed analysis, however, of the late-type giant 
population requires spectroscopy, owing to the large scat- 
ter in extinction near the Galacti c Center and the in trin- 
sic variations in giant star colors. iBlum "eTaTI (l2003h pio- 
neered work in this area, using spectroscopic and photo- 
metric observations of the most luminous giants and su- 
pergiants (50% complete at K$ ~10 mag.) to construct 
an H-R diagram for stars in the inner 5 pc. They report 
that the GC star formation rate in the central few parsecs 
is largely similar to that of the bulge. Specifically, they 
find that the majority of stars formed more than 5 Gyr 
ago, though they also find evidence that significant star 
formation also occurred during the past 100 Myr. How- 
ever, their conclusions are limited by their bright magni- 
tude limit, which samples only short-lived evolutionary 
stages, for which theoretical models are unc ertain. 

In t his paper, we build on the work of IBlum et al.l 
(2003) in an effort to better characterize the late-type 
giant population within the central parsec. We report 
deep photometric and spectroscopic observations of 329 
late- type giants in the GC, complete to 50% at Ks ~15.5 
mag. Our observations for the first time include the 
helium-burning red clump, as well as the red giant branch 
and asymptotic giant branch. This improved magnitude 
limit allows for the most robust picture of the Galactic 
Center star formation history to date, as these stars are 
much better understood th an the supe r giants and lumi- 
nous AGB stars studied bv lBlum et alJ ([2003). In §2, we 
present our observations and the techniques used to con- 
struct the Hertzsprung-Russell (H-R) diagram. Section 
3 discusses the resulting H-R diagram and implications 
for the GC star-formation history. We discuss our results 
in §4 and conclude with §5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spectroscopic Observations and T e ff 
Determination 

We observed eight Galactc Center fields between 
March and Septem ber 2006 using the integra l field spec- 
trometer SPIFFI (Eis enhauer et alJ l2003allbl ) in con- 
j unction with the MACAO adaptive optics module 
(jBonnet et al.ll2003l ) mounted on the SINFONI ESO VLT 
facility. We observed these fields for the dual purpose 
of identifiying main sequence B stars outside the cen- 
tral parsec (Martins et al. 2007b, in prep) and gathering 
spectroscopic observations for a large number of late- type 
giants. We chose fields with a radial distance from Sgr 
A* of less than 20 arcsec, selecting regions outside the 
minispiral to avoid nebular contamination in Brackett 7, 
and choosing fields north of Sgr A* to minimize sepa- 
ration from the AO guide star. We also avoided regions 
with extremely bright stars (K < 9). As stars this bright 
are relatively rare, the selection bias in avoiding these 
stars is negligible. If the fields had instead been chosen 
at random, we expect no more than a few K < 9 stars 
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Fig. 1. — NaCo Ks band mosaic used to derive photometry. The 
eight regions for which we have SINFONI spectra are overplottcd 
in black. The position of Sgr A* is shown with an *. 
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Fig. 2. — Histogram of stellar positions, corresponding to the 
stars detected in the spectroscopic fields marked in Figure [l] 



would be selected. 

We observed each selected SINFONI field in 50 x 100 
mas pixel mode, resulting in 0.2 arcsec resolution and a 
4.2" x 4.2" field of view. We simultaneously observed all 
fields in H and K band, leading to a spectral resolution of 
R k, 1500. Figure Q] displays the location of the observed 
fields. The total exposure time for all fields was 4200 
s, except for the fields at (17", 17"), (5", 11"), and (- 
8", 8"), which had total exposure times of 6600 s, 7200 
s, and 2400 s, respectively. The fields span a range in 
projected distance from Sgr A* of 4 arcsec to 26 arcsec. 
A histogram of the observed stellar positions is shown in 
Figure [21 

We reduced the raw data using a standard proce- 
dure to perform flat -fielding, sky subtracti on, and wave- 
length calibration (|Schreiber et al.l I2004D . Following 
extraction, we removed the stellar continua by divid- 
ing by a second degree polynomial fitted to the line- 
free regions of the stellar spectra. The normalized 
spectra allow us to directly compare the VLT spectra 
to the nor malized catalog spe c tra in iWallace &: Hinklg 
(|1997ft and iKleinmann fc Halll (|1986h . A typical resul- 
tant GC K band spectrum following normalization is 
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shown in Figure [3l Th ree template spectra from the 
IWallace fc Hinkla(|1997l ) spectral library are also shown 
for comparison. 

To classify the stars in our sample, we used the 12 CO 
2.2935 /im, v = 2 — rovibrational band head, the 
strongest feature in our spectra and a well-known T R ff in- 
dicat or, for a given luminosity class (|Kleinmann fe HaDI 
.1986). We use only this feature, as it is largely free of 
nebular and telluric lines, making it more reliable than 
other H— and K— band signatures in our spectra. We 
note that in contrast to the work of iBlum et all (ll996bL 
|2003| ). we do not need the H+K band H 2 feature to re- 
liably predict effective temperatures, as our sample con- 
sists only of giants with K$ >10.3 (see §2.2). Due to 
their rarity, stars sufficiently bright to be supergiant or 
long period variable candidates are not contained in our 
sample. 

To define a T e g - 12 CO index relation, w e computed the 
CO index defined bv lBlum et"afl (|1996bl ). Although not 
as well-known as other CO indices (e . g. equ iv alent width, 
photom etric indices of iFrogel et al.1 (|1978[ )), IBlum et "all 
(|1996bf ) showed that this index correlates well with other 
CO indices. Furthermore, it has the advantage of be- 
ing insensitive to small variations in the nearby contin- 
uum because it does not require a fit to the continuum. 
The index is defined as CO % = [(1 — F ba nd/ F cont ) x 
100)], where both bands are 0.015 /im wide and the con- 
tinuum and CO band centers are defined as 2.284 fj,m 
and 2.302 /xm, respectively. We note that although we 
adopt t he same definition for the CO index as Blu m etlrtl 
(199Gb), our spectra are normalized while theirs are not, 
and thus, our computed indices are not directly compa- 
rable to theirs. We calculated uncertainties in the com- 
puted indices, assuming the noise is dominated by photon 
statistics of the source and background and that the un- 
certainty in the band is approximately equal to that of 
the nearby continuum. 

Following IBlum "elaTl ((2003T) . we used CO strength to 
estimate effective temperatures for the cool giants in our 
GC sample. To define a CO index versus T e ff relation, we 
used a set of archi v al com paris on star spectra taken from 
IWallace fc Hmkiel (1199711 and IKleinmann fcHal (fj~986h 
with well-determined effective temperatures given in the 
literature. In an effort to avoid systematic errors in our 
derived index— T e g relation, we used a large number of 
references, assuming temperatures derived from a wide 
variety of methods. A summary of the comparison star 
data is given in Table [H and the resulting relationship 
between CO index and T e g i s shown in Figure [4] The 
IWallace fc Hmkiel (fl99l and IKleinmann fc^IalllpSfih 
relationships are in very good agreement with each other 
(for stars in common, ACO < 0.5%), suggesting that the 
CO index is largely independent of resolution or observ- 
ing system. The relation is tightest for CO < 17. As no 
stars in our GC sample have computed indicies exceding 
this value, we did not include comparison stars with CO 
> 17 in our fit. The resulting index-T e ff relation using 
both empirical data sets is: log T c ff = 3.7351 — 0.0060 x 
CO - 0.00040 x CO 2 , for all stars with measured indices 
17 > CO > 3. The errors in the derived coefficients are 
0.0092, 0.0022, and 1.2xl0~ 4 , respectively. 

In order to separate cool giants from warmer giants 
and main sequence stars, warm stars with CO < 3 are 
removed from the sample. Of the 355 detected stars in 



our spectroscopy, 329 were categorized as CO-absorbers, 
and we assigned effective temperatures to these stars us- 
ing the above index— T c g relation. Prior to temperature 
calculation, we cross-correlated all CO-absorbing stars 
with a CO-template in order to remove radial veloci- 
ties. Errors in temperatures were estimated based on 
the noise in each spectrum and the intrinsic dispersion 
in our index— T e g re lation (cr = 0.0050). As noted by 
iRamirez et~ai1 (|1997f ). errors determined in this way are 
strictly only lower limits. The spectral classification and 
analysis of the remaining early-type (non- CO-absorbing) 
stars i n our sample is less straight-forward (M artins et al.l 
2007a). and we defer analysis of these stars to a later pa- 
per (Martins et al. 2007b, in prep). 

2.2. Photometric Observations 

We made photometric observations using the imaging 
system NAOS/CONICA ( NaCo), consisting of the adap- 
tive optics system NA OS dRousset et al.l I2003T) and the 
NIR camera CONICA (Har tung et al.ll2003l ) at the 8.2- 
m UT4 (Yepun) of the ESO VLT. In April 2006, we col- 
lected several AO-corrected images in each of H and K$ 
band with a pixel size of 8x — 27 mas. The total ex- 
posure time was 64 s for each band, broken into thirty- 
two dithered images, in which every fourth image was 
separated from the previous exposure by 16 arcsec. We 
used our standard reduction pipeline to perform sky sub- 
traction, bad pixel corrections, flat field corrections, and 
stacking of images to create final mosaics. In Figure [TJ 
we display a sub-section of the final Ks band mosaic, 
containing all stars analyzed in this paper. Photometry 
was extracted for only this region to help ensure PSF 
constancy across the image. 

We used the crowded fie ld photometry package 
StarFinder (|Diolaiti et al.l2000l ) to establish relative pho- 
tometry and source detection. An empirical PSF for the 
central core (FWHM ~ 0.10 arcsec) was extracted from 
each image using seven bright, isolated stars. The full 
radial extent of the extracted PSF was r ~ 0.19 arc- 
sec. To derive the photometric curve of growth needed 
to place these results on an absolute scale , we a dopted 
the MTF-fitting technique of ISheehv et~"aTI (|2006l ). This 
technique fits the power spectrum of the image using a 
combination of the source spatial distribution function 
determined by StarFinder and a parameterized descrip- 
tion of the modulation transfer functions (MTF) of the 
atmosphere, telescope, AO system, and science camera. 
The advantage of this technique is that it derives the 
PSF encircled energy curve of growth, including the ex- 
tended, seeing-limited halo, and thus provides the aper- 
ture correction required for absolute photometric cali- 
bration. In deriving the curve of growth from the data 
themselves, we avoid systematic calibration errors intro- 
duced from using a PSF standard acquired at a different 
time and under different observing conditions than the 
target data. Given the rapid varia t ions in AO perfor - 
mancc ()Fitzgerald fc Graham! l200l IVacca et al.ll2007h . 
this technique is crucial for deriving accurate absolute 
photometry. 

The MTF-fitting technique has so far only been ap- 
plied to data obtained for the Keck Observator y LGSAO 
system (jSheehv et al.ll200& IVacca et al.l 12007! ) . To ap- 
ply this technique to the VLT images, we used the 
ISheehv et al.l (2006) software with appropriate input pa- 
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Fig. 3. — A typical stellar spectrum from the Galactic Center fields. The position of the 12 CO 2.2935 /im v=2- Q rovibrational bandhead 
and nearby continuum used to calculate the CO-index is shown as a solid line. Three template spectra from the Wallace & Hinklc (1997) 
spectral library are also shown for comparison. 

mirror actuator spacing. The software uses the IDL 
procedure MPFIT 7 to perform a Levenberg-Marquardt 
least-squares fit of the parameterized model MTF to the 
data. The best fits to the H and K$ band image power 
spectra are shown in Figure [5l The spatial frequency, i> n , 
is normalized relative to the telescope cutoff frequency in 
the image plane (D te ;/A). 

All images, even seeing-limited images, contain power 
up to the spatial frequency cutoff, which is D/X for a 
circular pupil. To faithfully record all spatial informa- 
tion allowed by the telescope aperture, the pixel sam- 
plin g frequency must t herefore satisfy 1/Sx > 2D/X 
(see lSheehv et al.1 (|2006| ) for a detailed discussion). The 
Ks band data are Nyquist-sampled: (\ / 8x) / (2D / X) 
= f. However, the H-band data are undersamplcd: 
(1 / 8x) I (2D I A) = 0.76. This aliasing causes power 
at frequencies 0.76 < v n < 1 to corrupt frequencies 
0.52 < v n < 0.76. The contamination of frequencies 
between 0.52 and ~0.7 is negligible, as the power aliased 
to these frequencies consists only of detector and back- 
ground noise. This power is orders of magnitude less 
than the non-aliased power at these spatial frequencies, 
and thus, has very little effect on the total power spec- 
trum. On the other hand, the aliased power for contami- 
nated spatial frequencies v n > 0.7 is of the same order-of- 




FlG. 4. — The 2.3 /im CO index (see text for definition) as a 
function of tempe r ature for t he comparison stars taken from the 
IWallace fc Hinkld (119971 ) and IKleinmann fc Halll |H986T) catalogs. 
The solid line is a least-squares fit to both data sets for stars with 
17 > CO > 3. 



rameters for the VLT telescope pupil size and geometry, 
the camera platescale and pixel size, and the deformable 



Available at |http://co w. physics. wise. edu^craigm/idl/ 
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TABLE 1 
Comparison Star Properties 



Name 


Sp Type 


Catalog 


CO 


ACO 




T e ff (K) 


AT cff 








(%) 


(%) 


(K) 


(K) 


HR 7001 

1111, I \J\J 1 


AO Va 


WH 


1.760 


0.004 


9420 


60 


HR 4295 


Al V 


WH 


1.434 


0.003 


9600 


25 


HR 5054 

lilt 


Al Vp 


WH 


1.823 


0.004 


8760 


260 


HR 6378 

ill b 1 U 


A2 V 


WH 


1.953 


0.004 


8690 


78 


HR 4534 


A3 V 


WH 


1.695 


0.005 


8857 


185 


HR 4357 

111 L TiOiJ I 


A4 V 


WH 


1.516 


0.007 


8243 


150 


HR 7557 

1 11 L 1 U'J 1 


A7 V 


WH 


1.147 


0.004 


7500 


200 


HR 2943 


F5 IV-V 


WH 


1.176 


0.014 


6532 


39 


HR 6927 

ill Lh \J\J^i 1 


F7 V 


WH 


1.219 


0.014 


6087 


22 


HR 4375 

111 L rtO 1 tj 


F8.5 V 


WH 


1.540 


0.015 


5676 


66 


HR 21 


F2 III 


WH 


1.563 


0.011 


6847 


137 


HR 403 

111L rrUO 


A 5 TTT-TV 

i\. 11 1 1 V 


WH 


1.828 


0.008 


8420 


360 


HR 1412 


A7 III 


WH 


1.554 


0.007 


7690 


320 


HR 1457 


K5+ III 


WH 


12.633 


0.023 


3866 


35 


HR 2985 

1111 U U 


G8 III 


WH 


5.540 


0.018 


5001 


56 


HR 3003 

1111 JUUO 


K4 III 


WH 


13.944 


0.023 


3961 


99 


HR 3323 


G5 Ilia 


WH 


3.045 


0.017 


5136 


88 






KH 


3.039 


0.007 






HR 4031 


F0 III 


WH 


2.145 


0.011 


6880 


150 


HR 4069 

1111 L ±\J\JiJ 


M0 III 


WH 


14.143 


0.021 


3730 


149 


HR 4517 


Ml III 


WH 


13.894 


0.019 


3828 


5.3 


HR 4883 

111L TOOiJ 


00 TTTn 

VJU lllL^ 


WH 


1.819 


0.014 


5747 


54 






KH 


1.377 


0.007 






HR 5017 

111 t, JU 1 i 


F3 III 


WH 


1.395 


0.017 


7141 


181 


HR 5340 

111b \Jtj^\J 


K1.5 IIIp 


WH 


10.997 


0.021 


4233 


55 


HR 6299 

111L KjAiOCJ 


K2 III 


WH 


8.653 


0.020 


4571 


100 






KH 


8.530 


0.009 






HR 6703 

111 L vf i Ui) 


OS 5 TTT 

VJ O . {J 111 


WH 


4.426 


0.020 


5011 


35 


HR 6705 

111 V, vJ 1 \j\J 


K5 III 


WH 


12.859 


0.023 


3934 


12 






KH 


12.569 


0.010 






HR 7fi35 

111 L 1 


MO— III 


WH 


13.035 


0.025 


3867 


50 






KH 


12.869 


0.010 






HR 7806 

111 L 1 UUU 


K2.5 III 


WH 


9.993 


0.022 


4286 


100 






KH 


9.785 


0.010 






HR 7886 


M6 III 


WH 


16.722 


0.025 


3243 


79 


HR 8317 


K0.5 III 


WH 


7.319 


0.024 


4658 


19 


HR 8694 


KO- III 


WH 


6.691 


0.021 


4830 


100 






KH 


6.463 


0.010 






HR 8905 


F8 III 


WH 


1.265 


0.014 


5942 


12 






KH 


0.797 


0.007 






SWVir 


M7 III 


WH 


18.350 


0.024 


2966 


36 






KH 


18.033 


0.008 






BKVir 


M7- III 


WH 


20.449 


0.022 


2944 


3d 






KH 


20.730 


0.006 






RXBoo 


M7.5-8 III 


WH 


18.262 


0.024 


2786 


46 
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magnitude as the non-aliased power, so the power spec- 
trum at these frequencies is heavily influenced by the 
aliasing. We, therefore, chose only to fit the power spec- 
trum for all spatial frequencies v n < 0.70. As only detec- 
tor and background noise contribute to the power spec- 
trum for spatial frequencies v n > 0.8, we are excluding 
very little information about the PSF. 

The relatively minor discrepency in Figure [5] between 
the model and data power spectra in K$ band is likely 
due to uncertainties in the deformable mirror influence 
function, as an explicit description of the VLT CILAS 
mirror influence function was not available. We, there- 
fore, used the approximat e influence fun c tion f or the 
Keck Xinetics mirror from Ivan Dam et al.l ((2004) . The 
overall power spectra fits are satisfactory, and we do not 
expect that the photometric accuracy of our technique is 
significantly affected by this choice. 

The Strehl ratios implied by the PSFs reconstructed 
from the fits in Figure[5]are 10% and 19%, for H and Ks 
band, respectively. To test the integrity of these fits, we 



divided the images into nine subimages and fit the power 
spectrum of each subimage. This process yielded fairly 
consistent results for all subimages, implying an error in 
the H and Ks band Strehl ratios of 0.5% and 1.0%, re- 
spectively. This agreement between subimages suggests 
that PSF variations due to variable distance from the 
guide star and variable exposure time in the final mosaic 
are only present at a low level. 

We determined the photometric zero poi nt from obser- 
vations of the near-IR standard star 9178 (jPersson et alJ 
Il998h . obtained on the same night. The measured at- 
mospheric extinction coefficients used to correct for the 
difference in airmass between our standard stars and sci- 
ence targets were 0.06 mag. airmass -1 for H and 0.07 
mag. airmass -1 for Kg- 

We derived photometric errors and completeness using 
the standard technique of inserting and recovering arti- 
ficial stars of specified magnitude. Completeness is not 
uniform, so we calculated errors and completeness sepa- 
rately for each spectral field. To avoid artificially crowd- 
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Fig. 5. — Comparison of the image po wer spectra to the b est-fit 
models based on the parameterization of Shcchy ct al. (2006). The 
H band data are aliased, and we therefore, fit only the power spec- 
trum for u n < 0.7 (see §2.2). The spatial frequency, u n , is normal- 
ized relative to the cutoff frequency in the image plane (Dt e ;/A). 
The Strehl ratios implied by these fits are 10±0.5% and 19±1.0% 
, for H and Kg band, respectively. 



K s photometric completion 
H photometric completion 
Spectroscopic completion 

Oix 

Ks 




10 12 14 16 18 

magnitude 

Fig. 6. — H and Ks band photometric errors and completeness 
based on the artificial star tests described in §2. The spectroscopic 
completeness is estimated based on the number of CO-stars for 
which we could extract spectra, compared to the total number of 
stars detected in the photometry. 

ing our images, we inserted only fifty stars at a time to 
each spectral subimage and repeated this procedure five 
times for each magnitude bin. We computed the pho- 
tometric error at each magnitude using a Gaussian fit 



to the difference between the input and recovered mag- 
nitudes. Figure [6] shows the average photometric errors 
and completeness at each magnitude bin. Also shown is 
the spectroscopic completeness, estimated by comparing 
the number of stars detected in the H and K $ band pho- 
tometry to the number detected in the S1NFONI data. 
The incompleteness in the SINFONI data at K s -13-15 
mag. is caused by crowding and confusion with nearby 
bright stars in the field of view. 



3. RESULTS 

3.1. Hertzsprung-Russell Diagram Construction 

Using the derived photometry and effective temper- 
atures presented in the preceding section, we are able 
to place our GC sample on the Hertzsprung-Russell 
(H-R) diagram. To convert the measured photometric 
magnitudes to luminosity, we use models that assume 
the K band magnitu des are measured in t he Johnson- 
Cousins-Glass system (iBessell fc Brettlfl988| ). Therefore, 
we first transformed the measured Ks magnitudes to 
the Johnson-Cousins-Glass s ystem, using the approxi- 
mate transformations given in ICarpenterl(|2001l) . assum- 
ing the Na Co and 2MASS JHKs colo r systems are iden- 
tical, as in lMcCaughrean et al.1 (|2004f) . We next adopted 
empirical bolometric corrections of iFluks et all (1994) 
for M stars and th eoretical corrections of Girardi 2005 
( |http://pleiadi/pd.astro.it| ) for hotter stars. For the the- 
oretical corrections, we assumed a solar metallicity and 
a surface gravity of log g = 2.0, though the results are 
largely independent of these assumptions. For all tem- 
peratures in the Galactic Center sample, we found sys- 
tematic differences < 0.05 mag. in BC^ for A[M/H] = 
1.0 dex and A logg = 1.0 dex, leading to systematic dif- 
ferences in the luminosity of Alog(L/L ) < 0.02. Ad- 
ditionally, we assumed a distance to the Galacti c Center 
of 8.0 kpc (|Reidlll993t lEisenhauer et al.ll2003d ). There 
has been much recent debate concerning the precise dis- 
tance to the Galactic Center. However, the choice of 
distance has only a small systematic effect on our de- 
rived luminosities; a change in distance of Ac? = 0.4 kpc 
induces a change in luminosity of Alog(L/L ) ~ 0.05. 
Finally, we corrected each star's luminosity individually 
for extinction using the photometric colo r, the derive d 
T e ff , and the interstellar extinction law of iRiekd (|1999h . 
derived from NICMOS observations of Galactic Center 
stars. All results are tabulated in Table [2] 



TABLE 2 

Properties of the Late- type Galactic Center Stars 



ID 


ARA 


ADec 


K 


fit 


H 




CO 


log T off 


<7logT e// 


A K 


<7A K 


log(L/L ) 


°"log(L/L ) 




(") 


(") 


(mag) 


(mag) 


(mag) 


(mag) 


(%) 




(mag) 


(mag) 




1 


6.26 


6.24 


15.275 


0.042 


17.419 


0.072 


8.28 


3.658 


0.006 


2.99 


0.15 


1.88 


0.06 


2 


6.32 


6.05 


15.259 


0.042 


17.316 


0.072 


8.30 


3.658 


0.007 


2.86 


0.15 


1.83 


0.07 


3 


5.67 


5.91 


14.898 


0.031 


17.051 


0.056 


6.93 


3.674 


0.004 


3.03 


0.11 


2.09 


0.05 


1 


5.67 


6.13 


14.239 


0.013 


16.344 


0.049 


8.61 


3.654 


0.006 


2.93 


0.09 


2.25 


0.04 


5 


5.51 


6.24 


15.304 


0.042 


17.438 


0.072 


7.01 


3.673 


0.005 


3.00 


0.15 


1.91 


0.06 


6 


5.02 


6.48 


15.531 


0.042 


17.826 


0.204 


7.15 


3.672 


0.008 


3.23 


0.37 


1.91 


0.15 


7 


4.86 


6.29 


15.131 


0.031 


17.274 


0.072 


7.01 


3.673 


0.007 


3.01 


0.14 


1.98 


0.06 


8 


4.99 


6.16 


15.327 


0.042 


17.406 


0.072 


7.34 


3.669 


0.005 


2.91 


0.15 


1.86 


0.06 


9 


4.75 


6.13 


15.503 


0.042 


17.678 


0.072 


8.76 


3.652 


0.006 


3.03 


0.15 


1.78 


0.06 


10 


4.35 


6.34 


15.474 


0.042 


17.135 


0.056 


6.82 


3.675 


0.004 


2.30 


0.12 


1.57 


0.05 


11 


4.35 


6.51 


15.240 


0.031 


17.387 


0.072 


6.69 


3.677 


0.004 


3.02 


0.14 


1.96 


0.06 



GC Top-Heavy IMF 

TABLE 2 — Continued 
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ID 



ARA 

(") 



ADcc 
(") 



K 

(mag) 



(mag) 



H 

(mag) 



(mag) 



CO 

(%) 



log T eff <r log Te 



A K 
(mag) 



<JA K 
(mag) 



log(L/L ) (Ti og ( L / L@ ) 



12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 



3.54 
3.56 
3.48 
3.08 
3.05 
2.81 
3.27 
2.78 
4.35 
4.78 
5.13 
4.97 
5.53 
5.10 
5.16 
5.45 
5.08 
5.13 
5.05 
4.56 
4.29 
4.59 
4.40 
4.24 
4.32 
3.00 
3.11 
3.62 
4.02 
3.62 
2.92 
-6.97 
-7.24 
-7.10 
-7.32 
-7.94 
-8.99 
-8.80 
-9.29 
-8.69 
-7.45 
-5.59 
-6.75 
-6.75 
-6.18 
-6.29 
-6.10 
-5.97 
-6.16 
-6.51 
-6.80 
-7.32 
-7.45 
-7.64 
-7.48 
-8.23 
-8.91 
-8.59 
-8.13 
-8.18 
-7.94 
-7.75 
-7.97 
-7.61 
-8.45 
-9.21 
-8.64 
2.94 
2.67 
2.54 
2.59 
3.21 
3.00 
5.62 
5.83 



6.56 
6.40 
6.13 
6.13 
5.83 
5.86 
5.62 
5.32 
5.18 
5.37 
5.51 
5.10 
4.54 
3.46 
3.67 
3.91 
4.18 
4.51 
4.70 
4.56 
4.29 
4.13 
3.48 
3.70 
3.02 
3.32 
3.83 
3.56 
3.97 
4.05 
4.13 
10.15 
10.31 
9.85 
9.88 
9.88 
10.12 
9.67 
9.37 
9.21 
9.37 
9.21 
8.80 
8.56 
8.18 
8.07 
7.86 
7.29 
6.43 
6.64 
6.86 
7.78 
8.29 
8.56 
8.99 
8.88 
8.64 
8.18 
8.13 
8.32 
8.45 
8.26 
7.51 
7.21 
6.99 
6.53 
7.34 
12.93 
12.74 
12.45 
12.04 
12.69 
12.64 
11.88 
11.42 



14.318 
13.803 
14.969 
12.132 
14.123 
15.867 
15.596 
15.115 
11.399 
14.438 
15.436 
15.341 
14.770 
14.078 
15.087 
15.316 
13.152 
13.238 
13.980 
13.382 
15.536 
16.137 
13.845 
13.676 
13.194 
12.562 
15.976 
15.106 
15.532 
15.160 
12.738 
15.030 
14.751 
15.112 
15.075 
15.494 
13.124 
15.374 
15.199 
15.156 
11.937 
10.319 
15.448 
15.948 
15.273 
15.706 
15.506 
15.037 
13.182 
15.628 
14.604 
13.190 
15.413 
15.467 
12.976 
15.498 
15.444 
12.292 
14.233 
15.497 
15.375 
15.616 
12.569 
14.774 
10.843 
12.543 
15.402 
15.314 
14.953 
15.140 
15.189 
16.003 
15.868 
15.501 
11.834 



0.017 
0.013 
0.031 
0.003 
0.013 
0.078 
0.042 
0.031 
0.002 
0.017 
0.042 
0.042 
0.031 
0.013 
0.031 
0.042 
0.005 
0.005 
0.013 
0.009 
0.042 
0.078 
0.013 
0.009 
0.005 
0.003 
0.078 
0.031 
0.042 
0.031 
0.003 
0.026 
0.026 
0.026 
0.026 
0.040 
0.005 
0.040 
0.026 
0.026 
0.003 
0.002 
0.040 
0.090 
0.040 
0.040 
0.040 
0.026 
0.005 
0.040 
0.022 
0.005 
0.040 
0.040 
0.005 
0.040 
0.040 
0.004 
0.010 
0.040 
0.040 
0.040 
0.004 
0.026 
0.002 
0.004 
0.040 
0.049 
0.027 
0.027 
0.027 
0.075 
0.075 
0.049 
0.003 



16.427 
15.707 
16.899 
14.305 
16.226 
17.802 
17.751 
17.234 
13.519 
16.735 
17.641 
17.335 
17.054 
15.666 
17.025 
17.369 
15.121 
15.345 
16.298 
15.542 
17.606 
18.223 
15.739 
15.619 
15.295 
14.834 
18.227 
17.234 
17.603 
17.431 
15.117 
16.988 
16.765 
17.048 
17.022 
17.295 
15.144 
17.408 
17.210 
17.093 
13.686 
12.222 
17.392 
18.090 
17.416 
17.964 
17.688 
17.188 
15.300 
17.772 
16.766 
15.195 
17.374 
17.477 
14.980 
17.511 
17.651 
14.300 
16.282 
17.501 
17.420 
17.814 
14.867 
16.954 
13.278 
14.852 
17.464 
16.938 
16.828 
16.947 
16.632 
17.874 
17.752 
17.546 
14.160 



0.049 
0.012 
0.056 
0.006 
0.042 
0.204 
0.204 
0.056 
0.003 
0.049 
0.072 
0.072 
0.056 
0.012 
0.056 
0.072 
0.006 
0.012 
0.049 
0.012 
0.072 
0.204 
0.012 
0.012 
0.012 
0.006 
0.204 
0.056 
0.072 
0.072 
0.006 
0.058 
0.058 
0.058 
0.058 
0.060 
0.008 
0.060 
0.058 
0.058 
0.002 
0.002 
0.060 
0.245 
0.060 
0.245 
0.060 
0.058 
0.009 
0.245 
0.058 
0.008 
0.060 
0.060 
0.008 
0.060 
0.060 
0.005 
0.047 
0.060 
0.060 
0.245 
0.008 
0.058 
0.002 
0.008 
0.060 
0.055 
0.055 
0.055 
0.046 
0.210 
0.210 
0.080 
0.003 



7.73 
9.40 
6.90 
11.62 
8.65 
8.77 
9.25 
9.20 
12.10 
9.71 
8.55 
10.09 
6.77 
8.31 
9.92 
10.07 
10.32 
10.91 
12.37 
8.63 
7.99 
5.27 
8.65 
7.61 
9.46 
3.70 
9.31 
5.89 
7.97 
9.04 
12.73 
8.98 
7.28 
5.96 
6.46 
4.70 
10.30 
7.61 
8.72 
6.25 
5.88 
11.08 
6.35 
6.47 
8.85 
8.91 
4.64 
4.32 
7.56 
9.21 
7.46 
8.64 
6.44 
6.86 
12.14 
5.75 
9.59 
8.80 
8.98 
7.83 
9.34 
6.29 
10.02 
5.67 
12.81 
10.86 
8.56 
5.66 
7.95 
8.00 
4.90 
7.56 
6.33 
9.55 
13.30 



3.665 
3.643 
3.675 
3.611 
3.653 
3.652 
3.645 
3.646 
3.604 
3.639 
3.655 
3.634 
3.676 
3.658 
3.636 
3.634 
3.631 
3.622 
3.600 
3.653 
3.662 
3.692 
3.653 
3.666 
3.642 
3.707 
3.645 
3.686 
3.662 
3.648 
3.594 
3.649 
3.670 
3.685 
3.680 
3.698 
3.631 
3.666 
3.652 
3.682 
3.686 
3.619 
3.681 
3.679 
3.651 
3.650 
3.699 
3.702 
3.667 
3.646 
3.668 
3.653 
3.680 
3.675 
3.603 
3.687 
3.641 
3.651 
3.649 
3.664 
3.644 
3.682 
3.635 
3.688 
3.593 
3.623 
3.654 
3.688 
3.662 
3.661 
3.696 
3.667 
3.681 
3.641 
3.585 



0.005 
0.005 
0.005 
0.005 
0.006 
0.005 
0.007 
0.007 
0.005 
0.004 
0.004 
0.005 
0.006 
0.005 
0.011 
0.007 
0.006 
0.005 
0.007 
0.004 
0.004 
0.003 
0.006 
0.005 
0.005 
0.004 
0.008 
0.004 
0.006 
0.005 
0.005 
0.008 
0.006 
0.006 
0.006 
0.006 
0.008 
0.009 
0.007 
0.006 
0.004 
0.008 
0.005 
0.008 
0.009 
0.007 
0.005 
0.008 
0.008 
0.008 
0.006 
0.006 
0.005 
0.006 
0.009 
0.008 
0.008 
0.005 
0.009 
0.007 
0.007 
0.006 
0.009 
0.007 
0.013 
0.008 
0.007 
0.004 
0.005 
0.006 
0.005 
0.007 
0.006 
0.006 
0.007 



2.95 
2.62 
2.70 
2.96 
2.93 
2.67 
2.99 
2.94 
2.87 
3.19 
3.08 
2.74 
3.22 
2.17 
2.66 
2.82 
2.69 
2.88 
3.16 
3.01 
2.89 
2.95 
2.62 
2.71 
2.91 
3.24 
3.13 
3.00 
2.89 
3.17 
3.23 
2.71 
2.82 
2.72 
2.73 
2.53 
2.77 
2.84 
2.79 
2.71 
2.44 
2.58 
2.72 
3.02 
2.98 
3.15 
3.10 
3.05 
2.97 
2.98 
3.03 
2.78 
2.75 
2.82 
2.70 
2.83 
3.06 
2.78 
2.84 
2.79 
2.83 
3.10 
3.19 
3.08 
3.32 
3.18 
2.87 
2.26 
2.60 
2.50 
2.00 
2.60 
2.64 
2.82 
3.14 



0.09 
0.03 
0.11 
0.02 
0.08 
0.39 
0.37 
0.11 
0.01 
0.09 
0.15 
0.15 
0.11 
0.03 
0.12 
0.15 
0.02 
0.02 
0.09 
0.03 
0.15 
0.39 
0.03 
0.03 
0.02 
0.01 
0.39 
0.11 
0.15 
0.14 
0.02 
0.11 
0.11 
0.11 
0.11 
0.13 
0.02 
0.13 
0.11 
0.11 
0.01 
0.02 
0.13 
0.46 
0.13 
0.44 
0.13 
0.11 
0.02 
0.44 
0.11 
0.02 
0.13 
0.13 
0.03 
0.13 
0.13 
0.01 
0.09 
0.13 
0.13 
0.44 
0.02 
0.11 
0.03 
0.02 
0.13 
0.13 
0.11 
0.11 
0.09 
0.40 
0.40 
0.17 
0.02 



2.26 
2.27 
1.93 
2.99 
2.30 
1.49 
1.71 
1.88 
3.23 
2.24 
1.84 
1.68 
2.22 
2.02 
1.76 
1.73 
2.53 
2.55 
2.30 
2.63 
1.74 
1.61 
2.28 
2.42 
2.63 
3.20 
1.61 
2.03 
1.74 
1.96 
2.82 
1.83 
2.05 
1.91 
1.91 
1.72 
2.57 
1.80 
1.81 
1.88 
3.07 
3.59 
1.76 
1.68 
1.85 
1.74 
1.94 
2.12 
2.73 
1.69 
2.19 
2.61 
1.79 
1.78 
2.53 
1.81 
1.79 
2.97 
2.21 
1.72 
1.73 
1.85 
2.97 
2.20 
3.60 
2.95 
1.76 
1.65 
1.86 
1.74 
1.62 
1.45 
1.56 
1.67 
3.11 



0.04 
0.02 
0.05 
0.01 
0.04 
0.16 
0.15 
0.05 
0.01 
0.04 
0.06 
0.06 
0.05 
0.02 
0.06 
0.06 
0.02 
0.02 
0.04 
0.02 
0.06 
0.16 
0.02 
0.02 
0.02 
0.01 
0.16 
0.05 
0.06 
0.06 
0.02 
0.05 
0.05 
0.05 
0.05 
0.06 
0.02 
0.06 
0.05 
0.05 
0.01 
0.02 
0.06 
0.19 
0.06 
0.18 
0.06 
0.05 
0.02 
0.18 
0.05 
0.02 
0.06 
0.06 
0.02 
0.06 
0.06 
0.02 
0.04 
0.06 
0.06 
0.18 
0.02 
0.05 
0.03 
0.02 
0.06 
0.06 
0.05 
0.05 
0.04 
0.16 
0.16 
0.07 
0.02 
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TABLE 2 — Continued 



ID ARA ADcc K a K H a H CO log T off cr logTe// A K o Ak 

(") (") (mag) (mag) (mag) (mag) (%) (mag) (mag) 



log(L/L ) (Ti og ( L / L@ ) 



87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 



5.94 
5.48 
5.18 
4.83 
4.78 
4.35 
3.70 
3.43 
2.38 
2.59 
3.16 
3.73 
4.51 
5.10 
5.67 
5.40 
5.16 
5.91 
5.08 
3.59 
3.29 
3.13 
3.08 
3.46 
2.94 
2.51 
2.67 
-4.86 
-5.40 
-6.18 
-6.59 
-6.78 
-7.48 
-7.70 
-7.07 
-6.86 
-6.56 
-6.21 
-5.16 
-4.89 
-5.18 
-5.83 
-6.56 
-7.21 
-7.51 
-7.75 
-8.23 
-8.29 
-7.80 
-7.70 
-7.32 
-7.07 
-6.83 
-7.02 
-7.26 
-5.99 
-5.89 
-5.45 
-5.26 
-4.86 
4.94 
4.18 
3.51 
5.94 
5.83 
5.59 
4.97 
4.32 
2.73 
2.84 
3.00 
3.35 
3.73 
4.67 
5.32 



11.12 
11.45 
11.53 
11.74 
11.53 
11.69 
11.42 
11.39 
10.96 
10.77 
10.69 
10.53 
10.85 
10.66 
10.80 
9.96 
9.88 
9.45 
10.18 
9.45 
9.40 
9.53 
9.83 
10.12 
10.07 
10.23 
9.53 
12.72 
12.39 
12.69 
12.77 
12.66 
12.26 
12.45 
11.93 
11.85 
11.99 
12.10 
12.20 
12.12 
11.85 
11.56 
11.45 
11.53 
11.74 
11.58 
11.45 
11.12 
11.18 
10.83 
11.29 
10.85 
10.88 
10.53 
10.50 
11.18 
10.66 
10.31 
10.18 
9.85 
16.85 
16.82 
16.69 
16.58 
16.28 
16.09 
16.42 
16.28 
16.60 
16.25 
15.98 
15.90 
15.96 
15.93 
15.77 



13.811 
13.479 
15.624 
15.531 
15.510 
14.863 
12.593 
13.861 
12.391 
13.164 
14.694 
12.842 
11.905 
13.593 
15.151 
15.114 
15.134 
15.730 
13.307 
13.281 
15.212 
15.129 
13.819 
13.302 
14.905 
14.381 
15.518 
15.647 
14.498 
15.679 
16.205 
16.404 
15.019 
15.224 
15.166 
15.556 
15.669 
15.638 
17.135 
15.325 
15.313 
14.756 
14.853 
13.517 
14.319 
15.081 
15.838 
15.889 
14.080 
15.377 
14.858 
15.836 
14.494 
14.002 
14.588 
11.861 
13.946 
15.461 
14.684 
12.872 
15.552 
14.313 
15.608 
15.878 
15.617 
15.837 
12.691 
16.030 
15.914 
13.465 
15.508 
15.638 
15.501 
15.194 
15.439 



0.011 
0.008 
0.049 
0.049 
0.049 
0.027 
0.003 
0.011 
0.003 
0.005 
0.019 
0.005 
0.003 
0.008 
0.027 
0.027 
0.027 
0.049 
0.008 
0.008 
0.027 
0.027 
0.011 
0.008 
0.027 
0.019 
0.049 
0.042 
0.019 
0.042 
0.068 
0.119 
0.033 
0.033 
0.033 
0.042 
0.042 
0.042 
0.131 
0.042 
0.042 
0.033 
0.033 
0.008 
0.019 
0.033 
0.068 
0.068 
0.011 
0.042 
0.033 
0.068 
0.019 
0.011 
0.019 
0.003 
0.011 
0.042 
0.019 
0.004 
0.033 
0.014 
0.033 
0.073 
0.033 
0.073 
0.003 
0.073 
0.073 
0.007 
0.033 
0.033 
0.033 
0.026 
0.033 



16.328 
16.289 
17.880 
17.641 
17.643 
16.784 
14.495 
15.757 
14.350 
14.799 
16.701 
14.355 
13.794 
15.708 
17.216 
17.121 
17.133 
17.717 
15.413 
15.303 
17.153 
16.961 
15.715 
15.255 
16.849 
15.823 
17.613 
17.899 
16.366 
17.406 
18.072 
18.180 
16.737 
16.988 
16.674 
17.281 
17.310 
17.358 
19.005 
17.257 
17.172 
16.508 
16.636 
15.305 
16.096 
16.952 
17.605 
17.677 
15.992 
17.235 
16.667 
17.186 
16.366 
16.085 
16.339 
13.400 
15.939 
17.386 
16.372 
14.806 
17.771 
16.510 
17.420 
18.001 
17.657 
17.933 
15.117 
18.088 
17.716 
15.255 
17.399 
17.629 
17.434 
17.162 
17.512 



0.046 
0.046 
0.210 
0.080 
0.080 
0.055 
0.006 
0.042 
0.006 
0.006 
0.046 
0.006 
0.003 
0.008 
0.055 
0.055 
0.055 
0.080 
0.008 
0.008 
0.055 
0.055 
0.008 
0.008 
0.055 
0.042 
0.080 
0.224 
0.048 
0.072 
0.224 
0.224 
0.048 
0.053 
0.048 
0.072 
0.072 
0.072 
0.252 
0.072 
0.053 
0.048 
0.048 
0.010 
0.045 
0.053 
0.072 
0.072 
0.045 
0.053 
0.048 
0.053 
0.048 
0.045 
0.048 
0.002 
0.045 
0.072 
0.048 
0.007 
0.193 
0.043 
0.058 
0.193 
0.058 
0.193 
0.005 
0.193 
0.058 
0.007 
0.058 
0.058 
0.058 
0.046 
0.058 



10.99 
11.75 
6.90 
7.38 
4.91 
8.03 
5.95 
8.13 
11.31 
6.53 
5.54 
8.50 
5.57 
9.17 
8.89 
6.77 
8.75 
6.20 
11.02 
11.49 
7.64 
6.68 
9.31 
14.05 
8.69 
5.62 
7.98 
7.68 
7.14 
6.76 
5.31 
6.30 
7.69 
6.94 
5.40 
6.11 
5.47 
6.14 
6.87 
5.98 
5.02 
7.64 
8.90 
9.99 
7.41 
8.21 
7.01 
6.84 
10.11 
6.32 
7.69 
5.97 
10.31 
9.61 
7.36 
10.47 
8.36 
5.11 
7.67 
10.37 
8.39 
9.01 
7.22 
5.31 
7.14 
6.11 
14.50 
6.94 
5.83 
8.70 
4.92 
7.04 
6.61 
6.63 
7.90 



3.621 
3.609 
3.675 
3.669 
3.696 
3.661 
3.685 
3.660 
3.616 
3.679 
3.690 
3.655 
3.689 
3.646 
3.650 
3.676 
3.652 
3.682 
3.620 
3.613 
3.666 
3.677 
3.644 
3.572 
3.653 
3.689 
3.662 
3.665 
3.672 
3.676 
3.692 
3.681 
3.665 
3.674 
3.691 
3.683 
3.690 
3.683 
3.675 
3.685 
3.695 
3.666 
3.650 
3.635 
3.669 
3.659 
3.673 
3.675 
3.634 
3.681 
3.665 
3.685 
3.631 
3.641 
3.669 
3.628 
3.657 
3.694 
3.666 
3.630 
3.657 
3.649 
3.671 
3.692 
3.672 
3.683 
3.564 
3.674 
3.686 
3.653 
3.696 
3.673 
3.678 
3.678 
3.663 



0.008 
0.007 
0.005 
0.005 
0.006 
0.005 
0.003 
0.005 
0.004 
0.004 
0.006 
0.005 
0.003 
0.007 
0.005 
0.006 
0.007 
0.009 
0.007 
0.006 
0.009 
0.005 
0.006 
0.007 
0.005 
0.005 
0.007 
0.007 
0.005 
0.005 
0.005 
0.008 
0.005 
0.005 
0.004 
0.004 
0.004 
0.006 
0.005 
0.005 
0.004 
0.004 
0.005 
0.005 
0.004 
0.007 
0.004 
0.004 
0.005 
0.005 
0.004 
0.006 
0.006 
0.005 
0.005 
0.005 
0.005 
0.004 
0.004 
0.005 
0.007 
0.006 
0.005 
0.005 
0.006 
0.005 
0.006 
0.008 
0.008 
0.006 
0.006 
0.005 
0.006 
0.004 
0.005 



3.49 
3.90 
3.18 
2.96 
3.02 
2.67 
2.67 
2.63 
2.66 
2.26 
2.83 
2.06 
2.65 
2.93 
2.86 
2.81 
2.77 
2.79 
2.88 
2.74 
2.70 
2.55 
2.61 
2.56 
2.69 
1.99 
2.92 
3.16 
2.60 
2.40 
2.62 
2.48 
2.37 
2.45 
2.09 
2.40 
2.29 
2.40 
2.61 
2.71 
2.62 
2.42 
2.45 
2.43 
2.46 
2.59 
2.46 
2.49 
2.62 
2.60 
2.51 
1.85 
2.55 
2.88 
2.43 
2.06 
2.77 
2.71 
2.33 
2.64 
3.10 
3.06 
2.52 
3.00 
2.86 
2.95 
3.24 
2.89 
2.52 
2.46 
2.66 
2.79 
2.71 
2.76 
2.89 



0.09 
0.08 
0.38 
0.17 
0.17 
0.11 
0.01 
0.08 
0.01 
0.01 
0.09 
0.01 
0.01 
0.02 
0.11 
0.11 
0.11 
0.17 
0.02 
0.02 
0.11 
0.11 
0.03 
0.03 
0.11 
0.08 
0.17 
0.41 
0.09 
0.15 
0.42 
0.45 
0.10 
0.11 
0.10 
0.15 
0.15 
0.15 
0.51 
0.15 
0.12 
0.10 
0.10 
0.02 
0.09 
0.11 
0.18 
0.18 
0.08 
0.12 
0.10 
0.15 
0.09 
0.08 
0.09 
0.01 
0.08 
0.15 
0.09 
0.02 
0.35 
0.08 
0.12 
0.37 
0.12 
0.37 
0.02 
0.37 
0.17 
0.02 
0.12 
0.12 
0.12 
0.09 
0.12 



2.56 
2.83 
1.86 
1.79 
1.90 
1.92 
2.90 
2.30 
2.78 
2.49 
2.13 
2.47 
3.18 
2.49 
1.85 
1.92 
1.83 
1.68 
2.51 
2.45 
1.81 
1.81 
2.27 
2.27 
1.89 
1.92 
1.76 
1.82 
2.07 
1.53 
1.45 
1.28 
1.75 
1.73 
1.65 
1.60 
1.53 
1.56 
1.03 
1.82 
1.82 
1.88 
1.80 
2.29 
2.08 
1.80 
1.48 
1.48 
2.14 
1.74 
1.87 
1.27 
1.94 
2.29 
1.96 
2.79 
2.31 
1.79 
1.87 
2.62 
1.80 
2.26 
1.59 
1.74 
1.72 
1.71 
2.76 
1.58 
1.51 
2.37 
1.76 
1.69 
1.73 
1.87 
1.78 



0.04 
0.04 
0.16 
0.07 
0.07 
0.05 
0.01 
0.03 
0.01 
0.01 
0.04 
0.01 
0.01 
0.02 
0.05 
0.05 
0.05 
0.07 
0.02 
0.02 
0.05 
0.05 
0.02 
0.02 
0.05 
0.04 
0.07 
0.16 
0.04 
0.06 
0.17 
0.19 
0.05 
0.05 
0.04 
0.06 
0.06 
0.06 
0.21 
0.06 
0.05 
0.05 
0.05 
0.02 
0.04 
0.05 
0.08 
0.08 
0.04 
0.05 
0.05 
0.07 
0.04 
0.04 
0.04 
0.01 
0.04 
0.06 
0.04 
0.01 
0.14 
0.04 
0.05 
0.15 
0.05 
0.15 
0.02 
0.15 
0.08 
0.02 
0.05 
0.05 
0.05 
0.04 
0.05 



GC Top-Heavy IMF 
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ID 



ARA 

(") 



ADcc 

(") 



K 

(mag) 



(mag) 



H 

(mag) 



(mag) 



CO 

(%) 



log T eff <r log Te 



A K 
(mag) 



<JA K 
(mag) 



log(L/L ) (Ti og ( L / L@ ) 



162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 
215 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
233 
234 
235 
236 



5.43 
5.78 
5.83 
5.99 
5.08 
4.99 
3.81 
3.38 
3.16 
2.92 
2.78 
2.54 
3.21 
3.24 
3.62 
4.10 
6.10 
5.75 
4.56 
4.00 
3.89 
3.67 
3.21 
3.00 
2.46 
4.75 
5.35 
5.21 
4.86 
4.62 
4.29 
4.13 
3.81 
3.75 
3.46 
4.62 
4.67 
5.48 
5.75 
5.83 
5.56 
5.32 
17.98 
17.42 
16.96 
16.63 
16.42 
16.42 
15.55 
15.47 
15.39 
15.15 
14.61 
14.61 
15.07 
15.61 
15.82 
16.09 
16.36 
16.96 
17.60 
17.12 
16.98 
16.52 
16.31 
15.74 
15.34 
15.42 
15.34 
15.93 
16.44 
16.98 
17.25 
17.39 
17.77 



15.88 
15.80 
15.52 
15.58 
15.50 
15.15 
15.31 
15.31 
15.15 
15.66 
15.34 
14.80 
14.42 
14.66 
14.80 
14.88 
15.07 
14.50 
14.23 
14.28 
13.77 
14.09 
13.47 
14.09 
13.45 
13.82 
13.53 
13.01 
13.04 
13.28 
13.42 
13.18 
12.96 
12.82 
12.18 
12.31 
11.96 
12.34 
11.96 
12.50 
12.80 
12.64 
18.79 
18.66 
18.36 
18.68 
18.63 
18.39 
18.36 
18.68 
19.06 
18.58 
18.50 
18.04 
18.09 
17.90 
17.93 
18.01 
17.82 
17.98 
18.09 
17.74 
17.58 
17.55 
17.28 
17.39 
17.39 
16.98 
16.71 
16.71 
17.15 
17.17 
17.04 
17.20 
16.69 



15.235 
15.236 
15.698 
16.165 
16.363 
15.308 
12.615 
15.383 
14.040 
15.052 
14.778 
13.166 
15.317 
15.590 
15.449 
15.905 
15.197 
10.803 
12.404 
15.637 
12.332 
15.154 
13.180 
15.071 
15.667 
15.838 
15.493 
13.400 
15.419 
15.571 
15.229 
15.177 
14.855 
14.958 
14.986 
14.028 
15.444 
14.066 
15.084 
15.007 
15.206 
15.004 
13.875 
12.761 
15.226 
14.970 
15.582 
15.578 
15.736 
16.367 
16.023 
14.538 
15.400 
12.841 
15.299 
15.055 
15.329 
15.249 
16.300 
13.942 
15.513 
14.996 
13.668 
15.812 
14.495 
15.275 
15.470 
16.262 
15.123 
13.527 
13.639 
14.927 
15.273 
14.481 
14.440 



0.026 
0.026 
0.033 
0.073 
0.076 
0.033 
0.003 
0.033 
0.008 
0.026 
0.026 
0.004 
0.033 
0.033 
0.033 
0.073 
0.026 
0.002 
0.003 
0.033 
0.003 
0.026 
0.004 
0.026 
0.033 
0.073 
0.033 
0.007 
0.033 
0.033 
0.026 
0.026 
0.026 
0.026 
0.026 
0.008 
0.033 
0.008 
0.026 
0.026 
0.026 
0.026 
0.008 
0.005 
0.021 
0.021 
0.027 
0.027 
0.027 
0.083 
0.061 
0.015 
0.027 
0.005 
0.027 
0.021 
0.027 
0.021 
0.083 
0.008 
0.027 
0.021 
0.005 
0.061 
0.015 
0.027 
0.027 
0.083 
0.021 
0.005 
0.005 
0.021 
0.027 
0.015 
0.015 



17.402 
17.399 
18.004 
18.239 
18.561 
17.408 
14.638 
17.319 
16.061 
16.976 
16.934 
14.934 
17.176 
17.545 
17.362 
18.090 
17.898 
13.327 
15.485 
17.484 
14.286 
17.022 
15.232 
17.010 
17.564 
17.896 
18.000 
15.575 
17.675 
17.613 
17.173 
17.135 
17.095 
16.816 
16.925 
16.714 
17.517 
16.048 
18.396 
16.824 
17.276 
16.793 
15.987 
14.750 
17.257 
16.876 
17.474 
17.470 
17.868 
18.105 
17.959 
16.357 
17.426 
14.809 
17.248 
16.851 
17.293 
17.087 
18.281 
15.841 
17.565 
16.761 
15.518 
18.124 
16.480 
17.219 
17.529 
18.066 
17.461 
15.422 
15.636 
16.440 
17.212 
16.331 
16.309 



0.058 
0.058 
0.193 
0.193 
0.216 
0.058 
0.005 
0.058 
0.043 
0.046 
0.046 
0.005 
0.046 
0.058 
0.058 
0.193 
0.193 
0.002 
0.007 
0.058 
0.005 
0.046 
0.005 
0.046 
0.058 
0.193 
0.193 
0.007 
0.058 
0.058 
0.046 
0.046 
0.046 
0.046 
0.046 
0.043 
0.058 
0.043 
0.216 
0.046 
0.058 
0.046 
0.041 
0.005 
0.058 
0.048 
0.058 
0.058 
0.202 
0.202 
0.202 
0.047 
0.058 
0.006 
0.048 
0.048 
0.058 
0.048 
0.204 
0.041 
0.058 
0.048 
0.007 
0.202 
0.047 
0.048 
0.058 
0.202 
0.058 
0.007 
0.007 
0.047 
0.048 
0.047 
0.047 



9.64 
8.57 
8.79 
6.62 
8.51 
10.85 
13.49 
8.35 
9.34 
8.92 
9.07 
5.00 
6.99 
7.46 
9.52 
7.66 
10.05 
16.04 
14.27 
4.80 
12.22 
9.60 
11.63 
7.33 
7.33 
7.32 
9.48 
11.59 
4.92 
7.31 
7.96 
5.54 
7.21 
6.62 
7.58 
11.82 
8.54 
9.55 
10.30 
7.83 
7.35 
7.93 
12.13 
7.67 
8.52 
8.24 
6.17 
8.39 
6.30 
5.70 
3.62 
7.22 
7.51 
5.63 
6.82 
4.29 
5.46 
7.03 
8.44 
9.37 
9.45 
8.00 
9.04 
5.87 
9.03 
6.77 
6.90 
6.45 
6.69 
10.90 
10.57 
7.26 
8.20 
8.79 
7.90 



3.640 
3.654 
3.651 
3.678 
3.655 
3.623 
3.581 
3.657 
3.644 
3.650 
3.648 
3.695 
3.674 
3.668 
3.642 
3.666 
3.634 
3.536 
3.568 
3.697 
3.602 
3.641 
3.611 
3.670 
3.670 
3.670 
3.642 
3.612 
3.696 
3.670 
3.662 
3.690 
3.671 
3.678 
3.667 
3.608 
3.655 
3.641 
3.631 
3.664 
3.669 
3.662 
3.603 
3.665 
3.655 
3.658 
3.683 
3.657 
3.681 
3.688 
3.708 
3.671 
3.667 
3.689 
3.675 
3.702 
3.690 
3.673 
3.656 
3.644 
3.643 
3.661 
3.648 
3.686 
3.648 
3.676 
3.675 
3.680 
3.677 
3.622 
3.627 
3.670 
3.659 
3.651 
3.663 



0.005 
0.007 
0.006 
0.008 
0.006 
0.007 
0.007 
0.006 
0.006 
0.008 
0.005 
0.005 
0.006 
0.006 
0.005 
0.005 
0.006 
0.007 
0.007 
0.003 
0.007 
0.005 
0.007 
0.005 
0.005 
0.004 
0.009 
0.006 
0.005 
0.005 
0.007 
0.006 
0.005 
0.006 
0.007 
0.007 
0.011 
0.009 
0.010 
0.007 
0.009 
0.006 
0.006 
0.004 
0.005 
0.005 
0.004 
0.007 
0.006 
0.004 
0.004 
0.004 
0.008 
0.003 
0.006 
0.003 
0.004 
0.004 
0.006 
0.004 
0.006 
0.004 
0.004 
0.004 
0.004 
0.006 
0.005 
0.007 
0.005 
0.006 
0.005 
0.004 
0.006 
0.005 
0.007 



3.00 
3.01 
3.22 
2.91 
3.07 
2.88 
2.68 
2.68 
2.79 
2.66 
3.00 
2.48 
2.59 
2.73 
2.63 
3.06 
3.78 
3.31 
4.22 
2.60 
2.62 
2.56 
2.78 
2.70 
2.64 
2.88 
3.51 
2.97 
3.20 
2.86 
2.70 
2.76 
3.15 
2.59 
2.70 
3.72 
2.88 
2.73 
4.68 
2.52 
2.90 
2.47 
2.86 
2.77 
2.82 
2.64 
2.65 
2.62 
3.00 
2.43 
2.74 
2.53 
2.83 
2.77 
2.73 
2.53 
2.76 
2.56 
2.75 
2.61 
2.83 
2.44 
2.54 
3.28 
2.74 
2.72 
2.89 
2.52 
3.30 
2.57 
2.73 
2.07 
2.69 
2.55 
2.59 



0.11 
0.11 
0.35 
0.37 
0.41 
0.12 
0.02 
0.12 
0.08 
0.09 
0.09 
0.01 
0.10 
0.12 
0.12 
0.37 
0.35 
0.02 
0.02 
0.12 
0.02 
0.09 
0.02 
0.09 
0.12 
0.37 
0.35 
0.02 
0.12 
0.12 
0.09 
0.09 
0.09 
0.09 
0.09 
0.08 
0.12 
0.08 
0.39 
0.09 
0.11 
0.09 
0.08 
0.01 
0.11 
0.09 
0.11 
0.11 
0.36 
0.39 
0.38 
0.09 
0.11 
0.01 
0.10 
0.09 
0.11 
0.09 
0.39 
0.07 
0.11 
0.09 
0.02 
0.38 
0.09 
0.10 
0.11 
0.39 
0.11 
0.02 
0.02 
0.09 
0.10 
0.09 
0.09 



1.85 
1.89 
1.78 
1.54 
1.46 
1.72 
2.61 
1.71 
2.25 
1.81 
2.05 
2.62 
1.74 
1.67 
1.61 
1.67 
2.16 
3.49 
3.28 
1.69 
2.75 
1.70 
2.50 
1.88 
1.61 
1.64 
1.95 
2.49 
2.01 
1.74 
1.79 
1.91 
2.14 
1.90 
1.90 
2.53 
1.76 
2.21 
2.55 
1.81 
1.90 
1.79 
2.23 
2.82 
1.82 
1.86 
1.69 
1.60 
1.76 
1.30 
1.62 
2.02 
1.79 
2.85 
1.81 
1.90 
1.86 
1.76 
1.36 
2.22 
1.67 
1.78 
2.31 
1.85 
2.06 
1.82 
1.80 
1.35 
2.12 
2.31 
2.34 
1.68 
1.76 
2.00 
2.06 



0.05 
0.05 
0.14 
0.15 
0.17 
0.05 
0.02 
0.05 
0.04 
0.05 
0.04 
0.01 
0.05 
0.05 
0.05 
0.15 
0.14 
0.02 
0.02 
0.05 
0.02 
0.04 
0.02 
0.04 
0.05 
0.15 
0.14 
0.02 
0.05 
0.05 
0.04 
0.04 
0.04 
0.04 
0.04 
0.04 
0.06 
0.04 
0.16 
0.04 
0.05 
0.04 
0.03 
0.01 
0.05 
0.04 
0.05 
0.05 
0.15 
0.16 
0.15 
0.04 
0.05 
0.01 
0.04 
0.04 
0.05 
0.04 
0.16 
0.03 
0.05 
0.04 
0.01 
0.15 
0.04 
0.04 
0.05 
0.16 
0.05 
0.02 
0.02 
0.04 
0.04 
0.04 
0.04 
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ID 



ARA 

(") 



ADcc 
(") 



K 

(mag) 



(mag) 



H 

(mag) 



OH 
(mag) 



CO 

(%) 



log T eff <r log Te 



A K 
(mag) 



<JA K 
(mag) 



log(L/L ) o- log(L/LQ) 



237 
238 
239 
240 
241 
242 
243 
244 
245 
246 
247 
248 
249 
250 
251 
252 
253 
254 
255 
256 
257 
258 
259 
260 
261 
262 
263 
264 
265 
266 
267 
268 
269 
270 
271 
272 
273 
274 
275 
276 
277 
278 
279 
280 
281 
282 
283 
284 
285 
286 
287 
288 
289 
290 
291 
292 
293 
294 
295 
296 
297 
298 
299 
300 
301 
302 
303 
304 
305 
306 
307 
308 
309 
310 
311 



17.31 
16.98 
17.52 
17.31 
16.39 
16.12 
16.04 
15.28 
15.31 
15.50 
15.12 
18.71 
18.58 
17.87 
16.98 
16.17 
15.93 
16.33 
16.60 
16.90 
17.09 
17.28 
17.55 
18.14 
18.31 
18.90 
18.50 
18.31 
18.71 
19.17 
17.87 
17.74 
17.15 
16.98 
16.69 
16.77 
16.39 
15.55 
15.42 
15.96 
15.47 
15.50 
15.85 
16.01 
15.82 
16.50 
16.74 
16.71 
16.71 
16.96 
17.17 
17.20 
17.58 
18.28 
18.33 
18.06 
17.77 
17.63 
-6.67 
-7.51 
-7.70 
-8.29 
-8.37 
-9.15 
-8.26 
-6.72 
-7.07 
-8.13 
-8.48 
-8.72 
-9.94 
-10.07 
-10.12 
-8.88 
-8.80 



16.69 
16.77 
16.07 
15.71 
15.63 
15.42 
15.77 
15.47 
16.09 
16.42 
16.36 
13.90 
13.80 
13.93 
13.61 
13.66 
13.45 
13.23 
13.45 
13.45 
13.18 
13.18 
13.31 
13.58 
13.50 
13.20 
12.58 
12.28 
12.26 
12.18 
12.64 
12.26 
12.18 
12.31 
12.18 
12.99 
12.61 
13.10 
12.50 
12.23 
11.77 
11.56 
11.34 
11.15 
11.10 
11.58 
11.58 
11.31 
10.99 
10.66 
11.20 
11.58 
11.29 
11.50 
11.69 
10.77 
10.75 
10.34 
18.85 
19.31 
19.12 
19.01 
18.82 
18.74 
18.23 
18.14 
17.77 
17.63 
17.74 
17.63 
17.90 
17.28 
17.50 
17.28 
17.15 



15.226 
16.383 
15.026 
13.073 
15.333 
14.285 
15.011 
15.066 
15.352 
15.382 
14.726 
14.621 
13.958 
14.910 
14.856 
16.296 
15.774 
14.934 
16.083 
15.023 
15.868 
14.923 
15.131 
15.066 
15.182 
11.054 
14.975 
12.522 
15.774 
12.940 
15.189 
13.417 
16.087 
15.043 
12.248 
14.951 
12.552 
15.324 
14.718 
13.624 
15.327 
15.479 
13.934 
15.087 
15.398 
15.047 
14.931 
15.040 
15.251 
16.019 
15.337 
15.334 
15.161 
15.082 
15.143 
14.916 
14.598 
15.301 
14.814 
15.600 
15.219 
14.783 
15.067 
11.437 
14.468 
13.651 
15.142 
15.598 
15.391 
15.467 
15.508 
15.566 
15.794 
15.160 
15.482 



0.021 
0.083 
0.021 
0.005 
0.027 
0.015 
0.021 
0.021 
0.027 
0.027 
0.015 
0.016 
0.009 
0.027 
0.027 
0.100 
0.064 
0.027 
0.064 
0.027 
0.064 
0.027 
0.027 
0.027 
0.027 
0.002 
0.027 
0.003 
0.064 
0.005 
0.027 
0.008 
0.064 
0.027 
0.003 
0.027 
0.003 
0.039 
0.016 
0.008 
0.039 
0.039 
0.009 
0.027 
0.039 
0.027 
0.027 
0.027 
0.039 
0.064 
0.039 
0.039 
0.027 
0.027 
0.027 
0.027 
0.016 
0.039 
0.014 
0.029 
0.014 
0.014 
0.014 
0.002 
0.011 
0.005 
0.014 
0.029 
0.029 
0.029 
0.029 
0.029 
0.060 
0.014 
0.029 



16.971 
18.381 
16.916 
15.095 
17.192 
16.163 
16.470 
16.950 
17.243 
17.339 
16.606 
17.520 
15.403 
16.616 
17.116 
18.025 
17.581 
16.765 
18.269 
16.766 
17.650 
16.713 
16.904 
16.834 
16.984 
12.618 
16.649 
14.414 
17.527 
14.683 
17.110 
15.210 
18.519 
16.845 
14.120 
16.704 
14.315 
17.308 
16.613 
15.550 
17.333 
17.545 
15.885 
17.174 
17.480 
16.877 
16.922 
17.024 
17.147 
17.658 
17.459 
17.229 
17.013 
16.835 
16.849 
16.814 
16.558 
17.319 
16.883 
17.609 
17.213 
16.753 
16.961 
13.493 
16.434 
15.524 
17.064 
17.594 
17.322 
17.360 
17.631 
17.617 
17.879 
17.095 
17.506 



0.048 
0.204 
0.048 
0.006 
0.048 
0.041 
0.047 
0.048 
0.048 
0.058 
0.047 
0.065 
0.009 
0.044 
0.055 
0.198 
0.065 
0.055 
0.233 
0.055 
0.065 
0.044 
0.055 
0.055 
0.055 
0.002 
0.044 
0.005 
0.065 
0.005 
0.055 
0.006 
0.233 
0.055 
0.003 
0.044 
0.005 
0.065 
0.044 
0.009 
0.065 
0.065 
0.045 
0.055 
0.065 
0.055 
0.055 
0.055 
0.055 
0.065 
0.065 
0.055 
0.055 
0.055 
0.055 
0.055 
0.044 
0.065 
0.042 
0.051 
0.042 
0.042 
0.042 
0.002 
0.044 
0.007 
0.042 
0.051 
0.051 
0.051 
0.051 
0.051 
0.196 
0.042 
0.051 



4.79 
6.51 
10.57 
10.61 
6.66 
8.51 
4.92 
6.97 
6.55 
7.03 
10.30 
11.38 
6.33 
4.50 
8.54 
5.75 
5.53 
6.30 
4.92 
6.39 
5.28 
6.56 
7.08 
7.71 
7.14 
11.69 
7.40 
10.74 
5.80 
11.04 
7.94 
9.23 
6.83 
8.59 
12.93 
6.96 
4.08 
4.81 
6.34 
8.26 
4.94 
6.01 
7.63 
9.49 
6.14 
8.00 
8.07 
9.80 
10.42 
5.89 
8.96 
5.43 
4.87 
6.09 
5.51 
10.69 
11.05 
8.23 
10.09 
4.45 
5.71 
8.22 
5.94 
8.46 
9.28 
9.17 
7.86 
5.03 
7.87 
7.25 
10.16 
7.43 
5.97 
8.64 
9.00 



3.697 
3.679 
3.627 
3.626 
3.677 
3.655 
3.696 
3.674 
3.679 
3.673 
3.631 
3.615 
3.681 
3.700 
3.655 
3.687 
3.690 
3.681 
3.696 
3.680 
3.692 
3.678 
3.672 
3.665 
3.672 
3.610 
3.669 
3.624 
3.687 
3.620 
3.662 
3.646 
3.675 
3.654 
3.591 
3.674 
3.704 
3.697 
3.681 
3.658 
3.696 
3.684 
3.666 
3.642 
3.683 
3.661 
3.661 
3.638 
3.629 
3.686 
3.649 
3.691 
3.696 
3.684 
3.690 
3.625 
3.620 
3.659 
3.634 
3.700 
3.688 
3.659 
3.685 
3.656 
3.645 
3.646 
3.663 
3.695 
3.663 
3.670 
3.633 
3.668 
3.685 
3.653 
3.649 



0.004 
0.008 
0.009 
0.005 
0.005 
0.008 
0.004 
0.005 
0.007 
0.004 
0.006 
0.005 
0.004 
0.004 
0.004 
0.006 
0.007 
0.004 
0.003 
0.004 
0.006 
0.004 
0.005 
0.004 
0.004 
0.005 
0.005 
0.004 
0.005 
0.007 
0.005 
0.004 
0.006 
0.003 
0.006 
0.003 
0.002 
0.003 
0.003 
0.004 
0.005 
0.008 
0.006 
0.008 
0.005 
0.010 
0.009 
0.005 
0.009 
0.006 
0.007 
0.005 
0.002 
0.005 
0.004 
0.008 
0.005 
0.007 
0.007 
0.004 
0.004 
0.005 
0.003 
0.003 
0.004 
0.004 
0.007 
0.004 
0.005 
0.004 
0.006 
0.006 
0.007 
0.005 
0.007 



2.45 
2.80 
2.57 
2.77 
2.59 
2.59 
2.02 
2.63 
2.64 
2.73 
2.56 
4.04 
1.99 
2.39 
3.16 
2.41 
2.53 
2.56 
3.10 
2.43 
2.50 
2.49 
2.46 
2.45 
2.50 
2.06 
2.31 
2.57 
2.45 
2.34 
2.67 
2.46 
3.44 
2.48 
2.48 
2.43 
2.48 
2.80 
2.65 
2.67 
2.83 
2.91 
2.72 
2.89 
2.93 
2.53 
2.77 
2.73 
2.58 
2.28 
2.95 
2.66 
2.61 
2.44 
2.38 
2.58 
2.66 
2.81 
2.85 
2.84 
2.81 
2.74 
2.66 
2.86 
2.71 
2.58 
2.67 
2.82 
2.68 
2.64 
2.93 
2.87 
2.94 
2.68 
2.80 



0.09 
0.39 
0.09 
0.02 
0.10 
0.08 
0.09 
0.09 
0.10 
0.11 
0.09 
0.12 
0.02 
0.09 
0.11 
0.39 
0.16 
0.11 
0.43 
0.11 
0.16 
0.09 
0.11 
0.11 
0.11 
0.01 
0.09 
0.01 
0.16 
0.02 
0.11 
0.02 
0.43 
0.11 
0.02 
0.09 
0.01 
0.13 
0.08 
0.02 
0.13 
0.14 
0.08 
0.11 
0.13 
0.11 
0.11 
0.11 
0.12 
0.16 
0.14 
0.12 
0.11 
0.11 
0.11 
0.11 
0.08 
0.14 
0.08 
0.11 
0.08 
0.08 
0.08 
0.01 
0.08 
0.02 
0.08 
0.11 
0.11 
0.11 
0.11 
0.11 
0.36 
0.08 
0.11 



1.79 
1.42 
1.72 
2.58 
1.75 
2.10 
1.70 
1.86 
1.76 
1.77 
1.85 
2.44 
2.07 
1.90 
2.10 
1.32 
1.58 
1.90 
1.70 
1.81 
1.54 
1.88 
1.76 
1.76 
1.76 
3.06 
1.75 
2.72 
1.54 
2.45 
1.79 
2.37 
1.78 
1.75 
2.70 
1.83 
2.89 
1.89 
2.03 
2.41 
1.90 
1.84 
2.32 
1.86 
1.87 
1.79 
1.93 
1.81 
1.64 
1.37 
1.81 
1.81 
1.88 
1.81 
1.78 
1.76 
1.91 
1.79 
1.93 
1.81 
1.91 
1.97 
1.90 
3.35 
2.05 
2.33 
1.82 
1.78 
1.72 
1.69 
1.69 
1.74 
1.72 
1.78 
1.69 



0.04 
0.16 
0.05 
0.02 
0.04 
0.04 
0.04 
0.04 
0.05 
0.05 
0.04 
0.05 
0.02 
0.04 
0.05 
0.16 
0.07 
0.05 
0.17 
0.05 
0.07 
0.04 
0.05 
0.05 
0.05 
0.01 
0.04 
0.01 
0.07 
0.02 
0.05 
0.01 
0.17 
0.05 
0.02 
0.04 
0.01 
0.06 
0.03 
0.01 
0.06 
0.06 
0.04 
0.05 
0.06 
0.05 
0.05 
0.05 
0.06 
0.07 
0.06 
0.05 
0.05 
0.05 
0.05 
0.05 
0.04 
0.06 
0.04 
0.05 
0.03 
0.03 
0.03 
0.01 
0.03 
0.01 
0.04 
0.05 
0.05 
0.05 
0.05 
0.05 
0.15 
0.03 
0.05 
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i n 


ADA 


ZAUCC 


1 V 




LJ 
11 






lf) g leff 


<rio g T eff 


Ak 




logO/ L Q}) 


°\ os (L/L 




(") 


(") 


(mag) 


(mag) 


(mag) 


(mag) 


(%) 




(mag) 


(mag) 






312 


-8 56 


17 31 


15.606 


0.029 


17.457 


0.051 


8.29 


3.658 


004 


2.56 


0.11 


1.57 


0.05 


313 


-8.05 


16.98 


15.477 


0.029 


17.505 


0.051 


7.95 


3.662 


0.006 


2.83 


0.11 


1.74 


0.05 


314 


-6.94 


17.12 


15.234 


0.014 


17.284 


0.051 


9.04 


3.648 


0.007 


2.84 


0.10 


1.80 


0.04 


315 


-6.59 


16.98 


15.524 


0.029 


17.516 


0.051 


5.68 


3.688 


0.005 


2.80 


0.11 


1.79 


0.05 


316 


-6.86 


16.71 


15.647 
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3.2. Observed Hertzsprung-Russell Diagram 

The GC H-R diagram is shown in Figure El Our data 
clearly show the red clump at log(L/i ) ~ 1.7, as well 
as the upper red giant branch / early asymptotic giant 
branch. There is an indication of the AGB bump at 
log(L/L Q ) ~ 2.3, and we also detect some lower red gi- 
ant branch stars at lower luminosities, though the ob- 
servations are highly incomplete in this region. The 
data set in Figure [7] allows for a more robust analy- 
sis of the Galactic Center star formation history than 
any other previously published data set. Photomet- 
ric st udies are limited to the modeling of luminosities 
alone <lFliekdll987t iNaravanan et al.lll996l:lDavidge et al.l 
1997t iPhilipp et all 119991 : lAlexander k. Sternberg! 119991 : 
Figer et al. 120041) . due to intrinsic variations in late-type 
giant colors and the large variation in Galactic Center ex- 
tinction. This point is illustrated in FigureEJ which com- 
pares the observed H-R diagram to a color-magnitude 
diagram (CMD) for the same GC stars. The RGB/ AGB 
and RC populations are more clearly distinguished in the 
H-R diagram than in the CMD, and there is less scatter 
in T e ff than H — K , due to variations in GC extinction. 
As a result, only the K-band luminosity function can 
reliably be modeled with broadband photometry alone. 
The observed H-R diagram in this study is also an im- 
provement with respect t o previous spectroscopic work 
(|Blum et al.lll996bl l2003h . due to our improved magni- 
tude limit (~5 mag. deeper), which allows for the anal- 
ysis of well-populated regions in the H-R diagram, in 
which the evolutionary models are fairly well understood. 

The errors in temperature and luminosity are shown 
with the data in grey in Figure [7] Average errors in 
temperature and luminosity for stars with luminosities 
1.6 < log(i/ Lq) < 2.0 (the red clump region) are shown 
at right. The temperature limit for the minimum CO 
index and the 50% completeness limit for the average 
extinction correction (Ax = 2.75 mag.) are also shown. 
Solar metallicity (in red, left) a nd metal-poor (Z=0 .008, 
in blue, right) isochrones from iGirardi et al.l (|2000D are 
overplotted. The isochrones show the wide range of ages 
spanned by the Galactic Center population. Figure[7]also 
demonstrates the age-metallicity degeneracy inherent in 
this part of the diagram. Stars in the same part of the H- 



R diagram may represent a metal-rich, intermediate-age 
(< a few Gyr) population or a metal-poor, old (> 5 Gyr) 
population. For stars younger than ^5 Gyr, the Galactic 
Cente r stellar population i s known to be ap proximately 
solar (jRamirez et al.l 120001: ICarr et al J 120001) . For older 
stars, the Galactic Center metallicity distribution is not 
well-known. We will address this uncertainty later in the 
paper. 

3.3. Deriving the Star Formation History 

To investigate the star formation history implied by 
our sample, we begin with a qualitative discussion, fol- 
lowed by quantitative analysis in §3.4-3.5. We start by 
considering three models based on candidate star for- 
mation histories presented in the literature. The first 
scenario consists of an ancient burst of star formation 
7.5-8.5 Gyr ago, similar to the s ingle stellar population 
of the bulge. iGenzel et~a l. (2003) compared this scenario 
to the K band luminosity function (KLF) of the inner 
parsec. The second model consists of constant star for- 
mation between 10 Myr and 10 Gyr ago. This m odel is 
based on the best-fit model of iFiger et al.1 (|2004t ). who 
considered the Galactic Center KLF within 40 pc of the 
center. The third m odel corresponds t o the best-fit star 
formation history of B lum et"aLl (|2003l ). who fit the H-R 
diagram of asymptotic giant branch and cool supergiant 
stars within the central 5 pc, using four specified age 
bins. All models are summarized in Table [3] 

For each scenario, we generated model H-R diagrams 
using the synthetic colo r-magnitude diagram com puta- 
tion algorithm I AC- Star ([Aparicio fc Gallardl200l l. The 
algorithm uses a Monte Carlo approach to compute com- 
posite stellar populations on a star-by-star basis, using 
a specified set of evolutionary tracks. The code accom- 
modates several additional inputs, including the initial 
mass function, star formation rate function, and chemi- 
cal enrichment law. To compare the generated models to 
the data, we added Gaussian noise based on the average 
observed errors in luminosity and temperature. We also 
randomly removed stars from the models according to 
our estimated spectroscopic completeness (Fig [6]) . 

The models are shown against our Galactic Center 
sample in Figure [9l A qualitative comparison suggests 
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Fig, 7, — H-R di agram for the GC stars with solar-metallicity (in red, left) and metal-poor (Z=0.008, in blue, right) isochrones from 
Girardi ct al. (2000) overplotted. The magnitude limit for 50% completeness (see Figure [6]l and the temperature limit for which stars no 
longer show CO (see Figure|4j are overplotted as solid dark lines. Errors are shown in grey. The mean errors in temperature and luminosity 
for stars with luminosities 1.6 < log(L/L0) < 2.0 (the red clump region) are shown at right. 



that the ancient burst model is insufficient to fully de- 
scribe the observed data set. The conti nuous star forma- 
tion model and the iBlum et "all (|2003h model span the 
full range in temperature and luminosity of the observed 
data set but seem to have relatively too many stars at 
cool temperatures. It is unlikely this observed effect is 
caused by systematic errors or a selection bias. A system- 
atic underestimate of the GC CO indices would shift the 
pattern to hotter temperatures. However, the red edge 
of the Hertzsprung gap (log T e g ~ 3.70) is well matched 
by the models, making such a systematic effect unlikely. 
In addition, some cool RGB/AGB stars are present in 
the data, while a systematic shift to warmer tempera- 
tures would allow none. Finally, the brightness and CO 
strength of cool RGB/AGB stars compared to warmer, 
dimmer red clump stars means that cool RGB/AGB 
stars are relatively easy to detect and classify. There is 
no obvious way of selectively removing these stars from 
the sample. 

Motivated by the comparison in Figure [9j we consider 
three alternative scenarios that could potentially explain 



the relative paucity of stars at low temperatures: 

1. The first and simplest possibility is that the star 
formation rate was low at early times (> 5 Gyr). 
To test this hypothesis, we generated two models 
consisting of linear combinations of constant star 
formation in two specified age bins (Models 4 and 
5 in Table [3|). Model 4 assumes bins of 0.01-5 Gyr 
and 5-12 Gyr. Model 5 assumes bins of 0.01-7 Gyr 
and 7-12 Gyr. 

2. A second possibility is that the oldest stars are 
metal poor. The previously described models as- 
sume a solar metallicity for all times . This assump- 
tion i s motivated by the work of iRamrrez et al.l 
(2000), who found the stellar [Fe/H] abundance to 
be approximately solar for stars younger than ^5 
Gyr. For older stars, the Galactic Center metallic- 
ity distribution has not been determined. However, 
the bulge, which formed 7-12 Gyr a go, is known 
to ha ve a nearly solar distribution (| Sadler et al.l 
1996). An extremely metal poor ancient popula- 
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TABLE 3 

Summary of Star Formation History Models 
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Two-bin, solar 
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0.04 
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Two-bin, solar 
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a Poisson likelihood parameter: Xa = ^ Si m i — n i + n i hi ! n i I s the number of observed stars and m; is the number of model stars in bin i. 



Goodness of Fit Parameter: Percentage of Monte Carlo trials in which synthetic data sets composed of stars drawn from the best-fit models have 
larger X\ than the Xa implied by the observed data set. If the model is an accurate representation of the true star formation history, this percentage 
should be approximately 50%. c Average relative star formation rate for each age bin specified in the third column. The rates are normalized such 
that the total relative star formation rate is 1.^ Uncertainty in the relative star formation rate, derived by fitting the star formation history to a 
series of synthetic data sets consisting of stars drawn from the observed data set. 



3. A final way to explain the relative lack of stars 
at cool temperatures invokes a non-standard ini- 
tial mass function (IMF). In all above models, we 
assume a Salpeter IMF with an upper mass limit 
of 120 Mq and a lower mass limit of 0.7 Mq. As 
stars less massive than ~ O.7M have main se- 
quence lifetimes comparable to the age of the Uni- 
verse, we are not sensitive to stars below this mass 
limit. Recent results indicate that the current ini- 
tial mass function of the Galactic Center is flatter 
than a Salpeter func tion or has an unusally high 
lower mass cut-off (|Nayakshi n~fc Sunyaevl 2005; 
iPaumard et"aHl2006D . Such a mass function could 
explain the relative lack of low mass stars in the ob- 
served old stellar population. To test this hypothe- 
sis, we consider three models, all assuming contin- 
uous star formation at solar metallicity for 0.01-12 
Gyr. The first model, Model 8, assumes a flat sin- 
gle power-law slope (dN/dm = ro~ , 85 ), ch osen to 
match the results of IPaumard et al.1 (|2006f h Mod- 
els 9 and 10 both assume a standard Salpeter slope 
(dN/dm = m~ 2 - 35 ), but with lower mass limits of 
2.5 Mq and 1.5 Mq, respectively. 

Results stemming from these three hypotheses are tested 
in §3.5. 

3.4. Quantification of Fit 

To quantitatively compare the models described above 
to the da ta, we adopted the numerical techniques de- 
scribed in Dolphin (2002). We first binned the observed 
and model H-R diagrams in temperature and luminos- 
ity for stars above our 50% completeness threshold, us- 
ing uniform bins with size three times our average errors 
(Slog(L/L Q ) = 0.12, 6logT eff = 0.018). The choice 
of bins in this technique is somewhat subjective. How- 
ever, tests using different bins sizes showed that while 
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Fig. 8. — Comparison of the observed GC H-R diagram to a 
color-magntidue diagram for the same GC stars. The figure demon- 
strates the advantange of the spectroscopic temperature derivation. 
The RGB / AGB and RC populations are more clearly distinguished 
in the H-R diagram than in the CMD, and there is less scatter in 
T e ff than H — K, due to variations in GC extinction. 



tion would, therefore, represent a completely sepa- 
rate population from the bulge. Given the current 
observed differences between the bulge and the nu- 
cleus discussed in §1, we can not exclude this pos- 
sibility. We, therefore, generated two models to 
represent this scenario. The first model, Model 6, 
consists of a linear combination of constant star 
formation in two age bins with solar metallicity for 
0.01-5 Gyr and a metallicity of 0.008 for 5-12 Gyr. 
The second model, Model 7, is a simple closed box 
model assuming constant star formation for 0.01-12 
Gyr, a starting metallicity of 0.004 and an ending 
metallicity of solar. 
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Fig. 9. — Comparison of the observed data (red crosses) to several proposed models given in the literature (black dots). See the text 
and Table [3] for a description of the models. We have added noise to the models based on the average observed errors in luminosity and 
temperature (upper left). We also corrected for the incompleteness function of the data (Fig [6}. A qualitative comparison suggests that 
all literature models produce relatively too many cool stars to describe the observed data set. 



the binning scheme does change the fit quality, it does 
not significantly affect the derive d star formatio n rates. 
This finding is in ag reement with [Dolphin (|2002D . 

For models with two age bins, we used the Num eri- 
cal Recipes procedure AMOEBA (jPress et al.lll992f ) to 
search for the linear combination of models that mini- 
mized the Poisson maxi mum likelihood parameter: \\ — 
2 Ei "*i - nt + m In (|Dolphinl ll99^. Here, rij is the 
number of observed stars and is the number of model 
stars in bin i. For models with fixed relative star forma- 
tion rates, we scaled the model distribution to minimize 
x\- To estimate the errors in the derived star forma- 
tion history , we im plemented the technique described in 
iBlum et alJ (|2003t ). We built a set of 100 H-R diagrams 
consisting of a random sampling of 329 stars drawn from 
the observed H-R diagram, allowing each observed star to 
be selected any number of times. We then re-derived the 
star formation history for each H-R diagram and given 
model. The resulting standard deviation in the derived 
star formation rate was taken as the la uncertainty. 

We measured the fit quality through a second set of 
Monte Carlo simulations. In this set, we generated 
10,000 synthetic data sets drawn from the fitted models, 
selecting 329 stars randomly for each trial and re-deriving 
the star formation history for the selected subset. For a 
model that is an accurate representation of the true star 
formation history, X\ derived from the actual data set 
should be comparable to xi. derived from a typical syn- 



thetic data set. Therefore, to establish the goodness of 
fit, we calculated the percentage of trials, Pa, in which 
x\ was larger when fitting the synthetic data sets to the 
models than when fitting the true data sample to the 
models. If the model is a good representation of the true 
star formation history, P\ should be approximately 50%. 

3.5. Model Results 

The results of all model fits are listed in Table El 
Hess diagrams showing the difference between the ob- 
served data histogram and the best-fit model histograms 
are shown in Figure 1101 White indicates regions where 
the model produces too many stars relative to the ob- 
served data set, and black indicates regions where the 
model produces too few stars relative to the observed 
data set. Examination of the Pa values in Table [3] con- 
firms the qualitative result discussed in §3.3 that no lit- 
erature model (Models 1-3) is a likely description of the 
observed data set. This result is also reflected in the 
first three panels in Figure fTUl which show white diago- 
nal streaks, corresponding to the overdensity of cool stars 
in the models relative to the data. 

Table[3]also suggests that the early-history low star for- 
mation rate models (Models 4 and 5) are unlikely repre- 
sentations of the data, though Model 5 returns a consid- 
erably better fit, and we are unable to completely exclude 
this model. However, examination of the fourth and fifth 
panels in Figure [TU] show that both models systematically 
overpredict the number of cool stars, as in Models 1-3. 
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Fig. 10. — Hess diagrams showing the difference between the observed data histogram and the best-fit model histograms summarized in 
Table [3] White indicates regions where the model produces too many stars relative to the observed data set, and black indicates regions 
where the model produces too few stars relative to the observed data set. The first three panels correspond to models proposed in the 
literature and described in the text. All show white diagonal streaks, corresponding to the overdensity of cool stars in the models relative 
to the data. The early-history low star formation rate models (Models 4 and 5) suffer similar systematic deficiencies. The early-history 
low metallicity models (Models 6 and 7) predict a red clump / horizontal branch morphology that is too blue compared to the observed 
distribution. The high IMF mass limit models (Models 9 and 10) predict a hotter and more luminous red clump than is observed. The flat 
IMF model (Model 8) corresponds to the best fit and shows minimal systematic trends compared to the other models. 



The early-history low metallicity models (Models 6 and 
7) also show systematic deficiencies when compared to 
the data. Both models predict a red clump / horizon- 
tal branch morphology that is too blue compared to the 
observed distribution. Similar discrepencies between the 
models and data are present in the high IMF mass limit 
models (Models 9 and 10). Each predicts a hotter and 
more luminous red clump than is observed. However, 
while Models 9 and 10 predict no old, cool stars, the 
flat IMF model (Model 8) predicts few old, cool stars, 
in good agreement with the observations. Model 8 also 
gives a satisfactory goodness of fit (P\ = 40%) suggesting 
that this model is a likely representation of the observed 
distribution. 

Using the models listed in Table El we tested for pop- 
ulation differences as a function of distance from the 
cluster center. Such differences within the central 1- 
2 pc are expected due to mass segregation, and sev- 
eral broadband photo metry studies have found indica- 
tions of such an effect dPhilipp et al.lll999i : iGenzel et al.l 
[2001 ISchoedel et al.H2007ft . To test for this effect in our 
data, we separated the data into two sets containing the 
four innermost (4" — 15", 11" median) and outermost 
(13" — 26", 20" median) spectral regions shown in Fig- 
ure [1] We then refit the two best-fit models (Models 5 
and 8) to each subset, applying the average errors and 
completeness specific to that subset. This process re- 
turned fits for each subset that are consistent with those 
obtained for the entire data set. In particular, the subset 
P\ values were within 10% of that obtained for the entire 
data set, and the relative star formation rates returned 



by both subsets for Model 5 were within la of the rates 
returned by the entire data set. We therefore, find no 
significant evidence in our data for a population gradi- 
ent. This result is no t that surprising, as th e broadband 
photometry results of lSchoedel et al.l (|2007l ) suggest that 
variations in the cool, low-mass stellar population within 
the central parscc should be most evident 3-7 arcsec from 
Sgr A*. Our sample is outside of the region expected to 
exhibit the largest population differences (see Figure [TJ. 
In the future, the technique we describe could be ap- 
plied to regions closer to Sgr A*, although limited spec- 
troscopic completeness due to increased stellar density 
would complicate such a study. Still, a thorough spec- 
troscopic study of the red-clump and RGB / AGB popula- 
tions as a function of distance from Sgr A* could provide 
the first defini tive test of mas s segregation within the 
central parsec (|Alexanderll2005Tl . 

With the restriction to simple models, the top-heavy 
IMF model appears to be superior to the remaining 
three model families we consider (low star formation rate 
(SFR) at early times, low metallicity at early times, IMF 
with high lower mass limit). The low-SFR, low-Z, and 
high Mi ower families all show consistent systematic de- 
ficiencies when compared to the data. It is, therefore, 
unlikely one could slightly change the chosen parameters 
(i.e. precise metallicity, age bin cut-offs, mass ranges, 
etc.) within these families to produce a model that is an 
adequate description of the data. 



4. DISCUSSION 
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Of the models considered in Table [3[ Model 8 (contin- 
uous star formation at solar metallicity with a top-heavy 
IMF) fits the observations best, and our Monte Carlo 
tests show that it is a reasonable description of the data. 
Based on the goodness of fit, we cannot completely ex- 
clude Model 5. However, as discussed in the previous 
section, this model shows systematic deviations from the 
data; we, therefore, favor Model 8. We note that Model 8 
is probably not the only possible description of the data 
with this degree of likelihood. All models so far discussed 
assume a simplistic description of the GC star formation 
history, and there are likely more complicated scenarios 
that also adequately describe the data. For instance, an 
additional model with age bins identical to Model 5 and 
an IMF slope identical to Model 8 yields a reasonable 
fit (xl = 257.0, P x = 26.3%). However, we do not feel 
that increasing the number of free parameters is justified, 
given the limits, uncertainties, and size of our data set. 

In the context of our current knowledge about the 
Galactic Center, continuous star formation with a top- 
heavy IMF is reasonable. At present, there is substan- 
tial evidence for stars in the Galactic nucleus spanning 
a wide range of ages, based on results from broadband 
photometry, and several distinct age tracers, including 
supergiants, AGB stars, and OH/IR stars (see §1). The 
present distribution and kinematics of gas in the inner 
Galaxy is also consistent with continuous star formation 
(jMorris fc Serabvn|[l996T) . 

In addition, the most recent epoch of star formation 
in the Galactic Center is l ikely represented by a top- 
heavy initial mass fun ction (|Navakshin fe SunvaevH2005l ; 
iPaumard et al.l [20061 ). If this recent epoch of star for- 
mation is not anomolous and periodic bursts of GC star 
formation have occurred throughout the history of the 
Galaxy, there is no a priori reason to believe the GC IMF 
would change with time. The present data set appears to 
support this picture. We note that the mass traced by gi- 
ants in the Central Cluster is primarily set by the initial 
mass formed, rather than dynamical effects, since mass 
segregation is expected to only str ongly affect the central 
~0.01 pc of the Galactic Center (|Hopman fc Alexander! 
l2006HFreitag et al.ll2006t) . On the 1-2 pc scale studied in 
this paper, the efficiency of mass segregation is thought 
to be much lower, in agreement with observational 
photo metry results (|Genzel et alj l2003t iSchoedel et al.l 
12007ft . 

Though our findings are broadly consistent with sev- 
eral stellar population studies of the Galactic Center and 
the nucleus at large, there remains some disagreement. 
In particular, our results are s omew hat different from 
those presented in Bl um et "ail (|2003l ). who probed the 
central 5 pc and used supergiants and bright AGB stars 
to quantify the star formation history in a method very 
similar t o that employed he re. Using their observed H-R 
diagram. iBlum et a l. (2003) argued for variable star for- 
mation over a wide range of ages, with the majority of 
stars formed more than 5 Gyr ago at solar metallicity. 
We find a much smaller fraction of our sample is repre- 
sented by solar-metallicity low-mass stars with ages > 5 
Gyr. 

The present work and the work of [Blum "et~ai1 (HH) 
study different regions in the H-R diagram. Each regime 
has relative advantages and disadvantages for studying 
the GC star formation history, and we discuss each in 



turn: 

1. Field of View: Due to the rarity of supergiants 
and bright AGB stars in the Galactic Center, 
IBlum et al.l ([2003) probed all late-type stars above 
their magnitude limit for the entire Central Clus- 
ter (r < 2.5 pc). They, therefore, derived absolute 
star formation rates and were able to compare the 
cluster mass implied to estimates from dynamical 
studies. In contrast, our study probes a relatively 
small region on the sky (~0.2 pc 2 within r < 1.0 
pc). As such, our results are somewhat susceptible 
to population inhomogeneities in the Central Clus- 
ter, and we are able to derive only relative star 
formation rates. 

2. Stellar Crowding: Spectral extraction from our 
SINFONI fields is complicated by stellar crowding. 
While we took care to ensure that our derived CO 
indices are largely independent of exact pixel and 
background selection, contamination by neighbor- 
ing stellar spectra remains a source of uncertainty 
for many o f the d im stars. The stars studied by 
IBlum et a l. ( 2003), on the other hand, are well sep- 
arated from each other. In addition, the GC AGBs 
and supergiants are much brighter than neighbor- 
ing stars, so there is negligible uncertainty in the 
Blu m et "ail (|2003f ) spectra due to crowding. 

3. Spectral Classification: The spectr a l clas sification 
method presented by Blum "etaTI (|2003) is very 
similar to that presented here. Our depper study 
has the advantage that only giants are observed, 
and thus, a separate luminosity class determina- 
tion is not required. Additionally, t he stars we 
study are warmer than those studied bv lBlum et alJ 
(2003), and thus, our CO-T c ff relationship is 
tighter than theirs, due to reduced uncertainties 
in mo del atmosphere spectr a at warmer tempera- 
tures (jReid^llawlcy 2000). We further note that 
the derived temperatures for t he red supergiants in 
the work of IBlum e t al. (2003) are likely systemat- 
ically too cool, as iLevesque et alJ ()2005l ) recently 
showed that red supergiants are ~400 K warmer 
than previously thought. 

4. Adequacy of Stellar Evolution Models: The evo- 
lutionary models in the part of the H-R diagram 
probed by our study (red clump, red giant branch, 
early AGB) are much less unce rtain than the part 
of the H-R diagram studi ed bylBlum et all (l2003f) 
(supergiants, TP- AGBs). iGallart et alJ ([20051 ) re- 
view the adequacy of AGB stellar evolution mod- 
els for deriving star formation histories, concluding 
that the observed bright AGB populations are of- 
ten more sensitive to poorly known modeling pa- 
rameters than the star formation history. The in- 
put physics for the late AGB stages is not well 
determined (e.g. mass loss, convection, efficiency 
of the third dredge-up), and curre ntly, only the 
Pado va libraries include these stages (jGirardi et al.l 
2000) . While uncertainties in mass loss and convec- 
tion are also present for the red clump, RGB, and 
early AGB, the input physics for these stages is 
much better understood, and a number of stellar 
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evolution models including these phases have been 
compared and tested against observations. For this 
reason, we believe our star formation history fit is 
more reliable than that presented in iBlum et al.l 
(2003). Finally, our data set also has the benefit 
that it contains several evolutionary features (red 
clump, AGB bump) that are cleraly distinguishable 
from surrounding regions in the H-R diagram. The 
diagram morphology assures us that no large sys- 
tematic effects are present in the data and allows 
us to construct physically motivated models. The 
region of the H-R diagram studied by IBlum et al.l 
(2003}) has no such morphological features. 

While we believe th at our results are more robust 
than those presented in lBlum et al.l (|2003h . we also note 
that the finding s are not necessaril y in conflict. The re- 
gion probed by IBlum et "ail (|2003f ) extends to ^2.5 pc 
from the center, whereas we probe only the central par- 
sec. Furthermore, we note that there are some simi- 
larities in our findings. Both studies find evidence for 
star formation throughout the history of the Galaxy, 
and both studies suggest that purely solar metallicity 
models are needed to produce the observed data. Still, 
further work is needed to resolve the remaining dis- 
crepencies (i.e., variable versus continuous star forma- 
tion, Salpeter versus flat IMF). Specifically, a sample 
tracing a large region on the sky and containing thou- 
sands of stars would represent a significant step for- 
ward in this field. The planned FLAMIN GOS-2 GC 
Survey on Gemini (jEikenberrv et all l2006f) will obtain 
4000 late-type giant spectra out to one degree in Galacto- 
centric radius (~140 pc) and with a spectral resolution 
of R^20,000. It will, therefore, provide unprecedented 
information on the Galactic Center chemical enrichment 
and star formation history. 

5. SUMMARY AND CONCLUSION 



We observed 329 late-type giants 4" — 26" north of 
Sgr A* with the integral field spectrometer SINFONI on 
the VLT. Combining spectral classifications with NaCo 
photometry, we derived luminosities and effective tem- 
peratures for these stars. Due to the improved mag- 
nitude limit of our sample relative to previous work, 
our derived H-R diagram clearly shows the red clump, 
as well as the red giant branch and asymptotic giant 
branch. Using a maximum likelihood analysis, we com- 
pared the observed distribution to models representing 
ten possible star formation histories. The best-fit model 
corresponds to continuous star formation over the last 12 
Gyr with a top-heavy IMF. The similarity of this result 
to the IMF observed for the most recent epoch of star 
formation is intriguing and perhaps suggests a connec- 
tion between recent star formation and the stars formed 
throughout the history of the Galactic Center. The up- 
coming FLAMINGOS-2 GC Survey on Gemini will pro- 
vide important information needed to further understand 
this suggestive result. 
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